Collisions of protons with light nuclei shed new light on nucleon structure 
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The high rates of multi-parton interactions at the LHC can provide a unique opportunity 
to study the multi-parton structure of the hadron. To this purpose high energy coUisions of 
protons with nuclei are particularly suitable. The rates of multi-parton interactions depend 
in fact both on the partonic multiplicities and on the distributions of partons in transverse 
space, which produce different effects on the cross section in pA collisions, as a function of the 
atomic mass number A. Differently with respect to the case of multi-parton interactions in 
pp collisions, the possibility of changing the atomic mass number provides thus an additional 
handle to distinguish the diverse contributions. Some relevant features of double parton 
interactions in pD collisions have been discussed in a previous paper. In the present paper 
we show how the effects of double and triple correlation terms of the multi-parton structure 
can be disentangled, by comparing the rates of multiple parton interactions in collisions of 
protons with D, and ^He. 
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The experimental evidence[T][5][3] jl] j5][6] and the beginning of the operations at the Large 
Hadron Collider (LHC) have recently triggered a lot of attention to the problem of Multi- 
ple Parton Interactions (MPI) in high energy pp collisions. Several papers have been writ- 
ten on the topic in the last few months and four international workshops have been orga- 
nized on the theme[7] |5] P] |TU] |1 1] |12] . Issues discussed in the literature range from estimates 
of the contributions of MPI in various reaction channels of particular interest for the LHC 
physicsp][ll][ig[ig|I3|l8||l9||20][2l][22][2g|2i| to the effects on the global features of the 
inelastic event and of the underlying event [25; \Z^r>\ |27] [28| |29] [30] [3T] [32] |33] . the QCD evolu- 
tion of the double parton distributions [34j[35j [3^ and the general formulation of MPI within 
QCD[37](M]|S]|lD]iII]il2]il3](ll](15J[46j[47j[48j[49j[53. Somewhat less attention has been de- 
voted to the study of MPI in hadron-nucleus collisions, although all effects of MPI are sizably 
enhanced in that case as a consequence of the much larger parton flux[51][52][53]. In our opin- 
ion, a good reason to pay more attention to pA collisions in this context is that, when studied 
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jointly with pp, MPI in pA collisions can provide a unique handle to study some aspects of the 
multi-parton structure of the hadron[5l]. 

In spite of being directly related to the multi-parton distribution functions, MPI in pp colli- 
sions can in fact provide only a partial information on the multi-parton distributions. Because 
of the localization of the large momentum transfer processes, the incoming parton flux, and 
thus the multi-parton distribution functions, depend explicitly on the relative transverse dis- 
tances between the interacting parton pairs |55)|56). At the same time, correlations in the hadron 
structure will prevent expressing the multi-parton distributions as an uncorrelated product of 
one-body distribution functions j57|. The rates of MPI will therefore depend both on the typical 
relative transverse parton distances and also on the moments of the multi-parton distribution 
in multiplicity. In the case of MPI in pp collisions, the two features are unavoidably linked 
in the measured cross section[5^|58) and, as a consequence, only a partial information on the 
multi-parton distributions can be obtained by measuring MPI in pp collisions. 

On the other hand, MPI in pA collisions can provide a further handle for a deeper insight into 
the correlated multi-parton structure |54j|58||. In pA collisions the MPI cross section is a function 
of the multiplicity of the target nucleons. In the case of two or more target nucleons, the dimen- 
sional scale factor, characterizing each MPI event, is provided both by the hadronic parameters, 
radius and partonic correlation length, and by the nuclear size. The hadronic scale measured 
with the generalized parton distributions [59j is rather small as compared to the nuclear scale, 
which thus acquires a dominant role even in the case of light nuclei. When two or more target 
nucleons take part to the hard interaction, the contributions to the double parton scattering 
cross section in pA collisions depend therefore only weakly on the hadronic dimensions and, by 
studying MPI in pA collisions, one can thus single out the effects of the moments in multiplicity 
of the multi-parton distributions from the effects due to the correlations in the transverse parton 
coordinates. 

Studying MPI in pA is simpler in the case of light nuclei, where the binding is not very 
strong: in this case the structure of the nucleon is not much affected by the binding and the 
non-relativistic form of the wave function in the rest frame of the nucleus is appropriate. In 
such a case it's simple enough to construct the boost, which allows to move from the rest frame 
of the nucleus, where the wave function is given, to the hadron-nucleon rest frame, which is 
most suitable to describe the collision and where the relativistic expression of the nuclear wave 
function is compulsory [60] |58] . Some relevant features of double parton collisions in proton- 
Deuteron interactions have been discussed in a previous article |5H]. In the present paper we will 
extend the analysis to the case of double and triple parton collisions of protons with Deuterium, 

and '^He; the main results presented in |58] will be reviewed and completed. 

The paper is organized as follows: the processes are described in a covariant way, using the 
formalism of Feynman graphs supplemented by effective vertices for the non perturbative dy- 
namics. The resulting expression is then reduced to a form containing the fractional longitudinal 
momenta and the transverse coordinates. The use of the non-relativistic nuclear wave function 
in the relativistic process is justified by the same argument of our previous study of double par- 
ton interactions in pD collisions [58j. The only difference concerns the technicalities, which are 
heavier when considering the three-body dynamics of a non-relativistic nuclear bound state (the 
problem is discussion in detail in Appendix A). Concerning the double and triple parton distri- 
butions, to remain as general as possible, we have not introduced any explicit expression with 
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correlation parameters. Rather we have hmited our discussion to the actual relations between the 
observables (namely the MPI cross sections) and non perturbative quantities, characterizing the 
double and triple parton distributions, directly related to the correlated multi-parton structure: 
Namely the various overlap integrals, with a strong dependence on the partonic correlations in 
transverse space, and the two functions of fractional momenta (one for the double and one for 
the triple multi-parton distributions) representing the deviation of the parton population from 
an uncorrelated, i.e. Poissonian, distribution. A simple and fully explicit correlated Gaussian 
model, where all quantities are worked out in detail, is presented in Appendix B. 

We like to anticipate a feature, which we will meet more than once and has no analogy in 
the case of MPI in pp collisions. The spread of the momenta of the bound nucleons will allow 
to produce the same initial partonic configuration in different ways. MPI in pA collisions is thus 
characterized by quantum interferences between initial state configurations, which differ in the 
nuclear fractional momenta and in the transverse parton coordinates. 

Since our main interest is to recognize the most important features of MPI in collisions of 
protons with light nuclei, we introduce a drastic simplification in treating all the particles entering 
the game as spinless bosons. In addition other finer details are neglected from the beginning: 
the proton and neutron masses are considered equal, as well as the binding energies of and 
^He, so p and p ^He collisions will be considered as equal. 

In Section 2 the double parton cross sections are worked out, distinguishing the cases with a 
different number of spectator nucleons. In Section 3 the triple parton cross sections are worked 
out, distinguishing the cases with the same pocedure. In Section 4 we discuss the relations 
between the accessible experimental information, namely the different cross sections, and the 
unknown quantities most directly related to the multi-parton correlations. The main points 
examined and the results obtained are finally summarized in the last part of the paper. 

2. DOUBLE SCATTERING ON DEUTERON OR TRITIUM 
2.1. Only one bound nucleon interacts with large momentum transfer 

Some aspects of double parton scattering of protons with Deuterium have been already dis- 
cussed inf58]. The process will be reviewed here and the analysis will be extend to the case 
of double parton collisions of protons with and ^He. In a double parton scattering on a 
Deuteron one has two possibilities, either only one nucleon interacts with large transverse mo- 
mentum exchange or the interacting nucleons are two. Analogously, in Tritium one may have 
either one or two spectator nucleons. With minor adjustments, the case of Tritium can thus be 
reduced to the case of Deuteron. 

The analytical expression for the hard scattering, when one of the component nucleons inter- 
acts twice and there are one (Deuteron) or two (Tritium) spectators, is conveniently expressed 
through the discontinuity of the forward scattering amplitude (see Fig.l). We start with the 
discontinuity of the amplitude J-2 of double scattering between two free nucleons: 
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Here and later on, we use the following notation: is the effective vertex for one-parton 
emission, <j)i the effective vertex for two-parton emission, i = p when the nucleon is the projectile 
proton, i = 1,2 labels the bound nucleons of the target nucleons. Tj is the T-matrix for the parton 
scattering (the index i labels the corresponding bound nucleon). The final state momenta are 
Qi, q'^, while Qi = Qi + 4% is the overall four-momentum of the final state i. The final directions 
are embodied in the angle Oj, the factor 1/8 relates the invariant relative phase space with the 
solid angle. We use the fact that the momentum variables have large and small components, so 
the "plus components of L, Ij, F3 are large and the corresponding minus variables are small. On 
the contrary the minus components of Ni,ai,Fi are large and the plus components are small. 
The four-momenta of the produced particles can have both plus and minus large components. 
More explicitly: large means oc ^/s, small means oc l/\/s, the transverse variables are constant 
with respect to centre of mass energy s and the general attitude, as in[58], is to integrate over 
the small components. 

It is useful to introduce now the fractional plus or minus momenta in the following way: 



L is the four momentum of the free proton, N the four momentum of the other nucleon, D 
the four momentum of the Deuteron, T the four momentum of the Tritium. The T-matrix 
amplitudes are related to the partonic cross section by: 



Hadronization is not included. The cross section is obtained from the discontinuity of the 
forward amplitude (L = L' , N = N'), removing the overall four momentum conservation and 
dividing by the incoming flux 2s. 

If one of the colliding nucleons is bound, then we define the fractional momentum of the 
nucleon as Z = 2N-/D-, in the case of the Deuteron, and Z = 3N^/T^, in the case of Tritium. 
It is useful to introduce also the variables Xj, defined as Xi = 2aj_/-D_, for the Deuteron, and 
Xi = 3ai_/r_, for the Tritium. The fractional momenta of partons with respect to the parent 
nucleons are thus zi = xi/{2 — Z), Z2 = X2/Z, for the Deuteron and analogously for the Tritium. 
In the expression of the flux is substituted by and T_ respectively. 

If both interacting nucleons are free, one obviously has = m^, A'^^ = m^. In the case of a 
Deuteron with a spectator nucleon, the relevant discontinuity is 



xi = X2 = I2+/L+, zi = ai_/A^_, Z2 = a2^/N. 



\T{1 + a ^ qq')\'^ = [Sirf {l+a^)daQ / dQ = {87rfxz{L+N^)d&Q/dn 



(2) 
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Figure 1: Double parton scattering in p D interactions. Only a single target nucleons interact with large 
monientun exchange 



Disc ^(2;o) = j dNdN' Disc T2iL, L' , D, D')5{{D - N) - {D' - N')) 



6{{D-Ny-m')/{2Try (3) 



where the Deuteron effective vertex <I?_d is introduced and the condition A'^^ = m? is substituted 
by L>2 = D'2 ^ ^2 

In the case of a or ^He target with two spectator nucleons, the corresponding discontinuity 

is 

Disc i3(2,o,o) = j dNdN'dNsBisc F2{L, L' , D, D')5{T -N-T' + N') 
<^t{N,N3)<^t*{N',N3) 



6{Ni - m')5{{T -N- N^Y - m') /{2T:f (4) 



where the Tritium (or ^He) effective vertex $t is introduced and the mass-shell condition is: 
Following [58], we proceed by defining the amplitude for finding one or two partons in the 
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projectile when the remnant of the parent nucleon has mass M. The integrated variable is 
A_ = — h)- and M± is the transverse mass. 
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The one-parton and two-parton amplitudes in the bound nucleon are defined in the same 
way. The only difference is that in the case of the bound nucleon one needs to replace 
with rn^ + N'^ . The covariant amplitude for finding a nucleon in the Deuteron is defined in an 
analogous way: 
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with the definition vE-^) (iV_ )/Ar_ = ^TdCI-D - A^)_)/(-D - Ar)_. We have also 

As a function of Zi, the amplitude ^]:)(N-)/N- has a maximum for Zi = 1, which obviously 
implies also Z2 = 1. 

Finally we define the covariant amplitude for finding two nucleons in the Tritium 
Trip / [iV2 - m2] • [iV| - m2] • [7V| - m^] H '^^.?+^(^+ " E 



2(27ri 

9 <I>T 



2T2 M^ZiZsZa/S - mj^ -^^2^3 - 2Z3Z1 - ^gZiZs 
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The expression is evidently symmetrical in (1,2,3) and can thus be indenti- 
fyed with ■^t{Ni_, N2-)/iNi_N2-), or with '^TiN2-, N3_)/{N2-N3_), or with 
^t(-^3-) -^i-)/(-^3_A^i_). The corresponding Tritium amplitudes have also a maximum 
in Zi, where Zi = 3, close although not strictly equal to one. A more detailed inspection 
shows that the lack of strict equality is due to the possible difference in the three transverse 
momenta. By defining 7j = 3m'^j_/M!^ — 1/3, we find that, at the first order in 7^, the Tritium 
amplitude is maximized for Zi = 1 — 7j/6 + 7i/2. 

We proceed now with the integration on the transverse variables, in the frame where the 
external transverse momenta L±, and D±, are equal to zero. We take the two-dimensional 
Fourier transforms ( bi is conjugated to Oj, Bj to Nj, /3i to li, all the variables are two-dimensional 
vectors) . 

■01 = {2^)"^ J 'tjji exp[ilb]db 1^2 = {^tt)"^ J 'tij2{bi,b2) exp[ilibi + il2b2]dbidb2 
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* = i"^^)'^ J ^iB)exp[iNB]dB (8) 

and analogously for the complex conjugated functions, with the variables b'^,b'2,B'. 

The integration over the transverse-momentum variables gives the diagonal property bi = b'l 
and so on. Moreover one obtains the geometrical relation: 61 — 62 = /?! — /32- 

The one-body and two-body parton densities are defined by the following integrals on the 
invariant mass of the residual hadron fragments: 

r(xi,X2;6i,62) = ^,77^ / \^PM{xl,X2■,b,,b2)\^-^'^^^^LldM^ (9) 

2{2lT)° J 1 — Xl — X2 

The residual dependence on gi_L, q[j_ is transformed into an angular dependence on Qi, 0,2- 
As an effect of the nucleon motion, 6)^ is coupled to the interactions by the integration 

on the fractional momentum Z, while the integration on the transverse variable B is decoupled 
from the other transverse variables. So the cross section is readily expressed as 

<i =J^sj r(^i,^2;6i,62)^^^^ —r{x,/Z,X2/Z;P^,/32) 

\^d{2- Z;B)\y{2- Z)dBdbi db2 dPi dP2 S{bi - 62 - /3i + P2) 

dxidx2 dx'i dx'2 dZ dO, \ dQ.2 (10) 

Note that the main dependence on Zi, for Z^ close to one, is due to the nuclear wave function 
while the error, due to the approximation Z = \ m the parton distribution F, is less 
important. The expression above can thus be slightly transformed as follows. 
From the properties of ^ we have: 



\^>d{2 - Z)|V(2 -Z) = [l + {1- Z)]\^>d{Z)\^/Z'' 

the second addendum is odd for the substitution Z ^ [2 — Z) so that the integration in Z which 
runs from to 2 gives zero and the cross section is more conveniently expressed by: 
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\^d{Z; B)\yz^dB dbi db2 dPi d^ S{bi - 62 - /3i + P2) 
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2.2. Two bound nucleons interact with large momentum transfer 

In colhsions of protons with D or ^H/^He, the presence of the nuclear wave function induces 
the presence of two kinds of contributions: "diagonal terms" in direct correspondence with the 
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Figure 2: Double parton scattering in p D interactions. Both target nucleons interact with large momen- 
tun exchange 

processes taking place when the nucleons are free and a number of nondiagonal or "interference" 
terms, which are due to the presence of the nuclear wave function^. 

The simplest case, where diagonal and nondiagonal terms appear, is the double scattering 
on a Deuteron affecting both bound nucleons. In this case the "diagonal" discontinuity has the 
already given form: 



Disc^ 
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[1] The discussion in[58' is limited to the cases where the contribution of the interference term is not relevant. 
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whereas the interference term has the form: 
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The diagonal term was akeady elaborated in|58j, so we are interested in the differences be- 
tween the two cases. In the diagonal case the conservation of the large components of momenta 
implies that they are equal on the two sides of the diagram = a_ = a'_, = N'_; 
the transverse variables become diagonal through the Fourier transformation. In this way the 
whole expression of the cross section is expressible in terms of densities i.e. square of the partonic 
wave function and square of the nuclear wave function. The cross section is thus expressed again 
through the one-body and two-body partonic densities: T{z; b), r(xi, X2; (32), obtained from 
the effective vertices cp, (p. The nuclear density is simply given by 

In the interference case the conservation of the large components of momenta still implies 
the equality on the two sides of the diagram 1+ = I'j^, a_ = a'_ but for the nuclear variables 
(iV — 02)- = {N' — ai)-, moreover the transverse variable b±, conjugated to a±, does not become 
diagonal through the Fourier transformation. The interference term cannot be expressed only in 
terms of partonic densities, we need to introduce a more complicated expression: 

1 /" ^2 T, ^2 ^1^2 



WiiZ, Z';xi,X2;bub2,B) = j^ J dUfdM^ 



[Z - xi){2 - Z' - X2) 

x^Mi(Si/^;6i)V'k(^i/(2 - Z');bi - B)ijM,{x2/i2 - Z)-b2)rM,{x2/Z';b2 + B) (14) 

The previous relation for the nuclear variables can be written also Z — Z' = X2 — xi. It is 
possible to factor the expression W symmetrically into two parts Wi = H(Z,Z';xi,X2',bi,b2 + 
B) X H{2 - Z,2- Z'; X2, xi; 62, ^1 - B) with 

H{Z,Z';xi,X2;bi,b2) = :^ [ dMf^===^^===i,MA^i/Z-M)rMA^2/Z'-,b2). 

2{2-kY J ^ {Z - xi){2 - Z' - X2) 

The part of the cross section coming from the direct term is (as given in|58j) by 

pD\ 1 f J., a n-,da{xi,xi)da{x2,X2) ,_ 
'^2,2\d = J r(xi,X2;/3i,/32) r(xi/Z;6i)r(x2/(2-Z);62) 

X \^D{Z;B)\^dBdbidb2d(3id^2S{B-bi + b2-(3i+^2) 

X Z~'^ dxidx2dxidx2 dZ dili dil,2 (15) 
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Figure 3: p D interactions with two target nucleons involved. Configurations in transverse space of the 
diagonal term 

In fig. 3 we show the corresponding configuration in transverse space. 
The part of the cross section coming from the interference term is 



In fig. 4 we show the configurations in transverse space of the two interfering amplitudes. 

A comparison between the two expressions: As far as the longitudinal variables are concerned 
the interference term requires the nuclear wave function to be taken at different values of Z, so 
it is depressed with respect to the diagional term. It must be taken into account that the width 
of the nuclear wave function is determined by the binding energy while the mismatch of the Z 
terms depend on the difference in fractional momentum of the partons. At large total energies 
(beyond the TeV) the values of x can be small, still mantaining the process within the limits of 
perturbative dynamics, so the depression is not necessarily very strong. For what concerns the 
transverse variables the possible depression is given by the ratio between the hadronic size (as 
seen in the hard QCD processes) and the nuclear size. The first scale sets the size of the variables 
bi, Pi the second characterizes the size of B. After integrating over B with the (5— function, one 
has B = (bi — b2 + Pi — ^2)- A simplified form for the interference term is thus obtained when 
neglecting the hadronic size as compared with the nuclear size by setting B = in nuclear 
wave function ^. Note that the effect depends on the transverse degrees of freedom, while the 
dependence on the total energy is weak, once the region of perturbative dynamics is reached. A 
quantitative, but model dependent discussion of this feature can be found in Appendix B. 

Conversely the denominators Z, Z' could be set equal to 1 in the factors F since the range of 
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Figure 4: p D interactions with two target nucleons involved. Configurations in transverse space of the 
two interfering amplitudes. Both configurations generate the same partonic initial state. 

variation of Xi is large, as compared with the variation of Z, allowed by the nuclear function 
From the previous treatment we learned that the presence of a spectator has a modest influence 
on the process, so the double scattering on a Tritium has only minor differences in comparison 
with the scattering on a Deuteron. The property holds also in the next cases to be considered, 
with the exception of those processes, which are possible in presence on Tritium but not on a 
Deuteron. 



3. TRIPLE SCATTERING ON DEUTERON OR TRITIUM 
3.1. Only one bound nucleon interacts with large momentum transfer 

The analytical expression for the hard scattering, where one of the component nucleons in- 
teracts three times and there are one (Deuteron) or two (Tritium) spectators, is again related to 
the amplitude J^s for the triple hard scattering between two free nucleons: 

n- X- 1 / 0p (Pl 01 , '\rr*t,' I . ^ 

'''''' = WWJ fA^l a'la'laf'^^'''' ^^'""'^ ''''^ 

X T2{l2a2 q2q2)T2* {I' 2a 2 ^292)^3(^303 9393)^1* (^'303 QsQa) 

X 5{L-"^l- Fi)5{L -^l' - Fi)5{Ni - ^ a - Fi)5{N'i - ^ a' - Fi) 

X 5{li + ai- Qi)5{l'i + a' I - Qi)5{l2 + 02 - Q2)5{^'2 + a'2 - Q2) 

X 8{h + as - Qi)5{h + a'3 - Q3)5{Ni - N'l) 

X JJd(a/8) (fai(fa'idHid%(fFj5{Fj^ - M/) d^Nid'^QidMj^ (17) 
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Figure 5: Triple parton scattering in p T or p '^He interactions. Only a single target nucleons interacts 
with large monientun exchange 

Analogously to the case of double scattering, also for the calculations of J's we use the property 
that, in the regime of interest, the momentum variables have large and small components; the 
term (p represents the vertex for emission of three partons from a nucleon. 

The factor: l/(i^Z|^|) is thus integrated in li_ and I2- independently from the rest of the 
diagram, in fact in all other conservation relations, the small terms li_ and I2- enter together 
with large components e.g. of Oj and can thus be neglected. Implementing the conservation 
(L — Fi) = li +I2 + h, the integration gives 

1 f 1 „ 4) I 4> 



We use here also the fractional momenta: Xi = li^/L^ ~ Qij^/Lj^ and, through the conservation 
for the plus components, we obtain: 

V's.mC ) 2L+'^ xiX2h\+ X2X'ili\+ X'iXil2\ + [M\/ {I -Y,Xi) - m?]xi,X2,x^ ^ '* 

M_|_ is the transverse mass of the remnants: M_|_^ = Mp"^ + , F± = — X] Zj_L 

The factor (p/ (af a| a|) can be integrated in the same way. One needs only to keep into 
account that we are interested in a situation where A^i is a generic time-like four vector, with 
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positive energy, no longer subjected to the condition N'^ = rr? . Even when Ni enters in a 
wave function, it is still almost on shell, as we are dealing with a weakly bound systems (it has 
nevertheless a transverse momentum Ni±). Let F be the four-momentum of the remnants; then 
in strict analogy with the previous result we get: 



M± is again the transverse mass of the remnants, here the conservation of the minus component 
is used. The integration on the remnants F can be treated as in the double scattering case: 
/ d^F5{F'^ - M'^)dM'^ = J dF±/F±d'^F±. The longitudinal integration is then performed by 
means of the 5-functions with the results: 1/Fi_ = l/[iV-(l — ^ -z)] l/-^4+ = ^ S 2;)]. 

Looking at the cut graph in Fig.l, we see that the equality of N2 and, in case, of N3 forces 
also Ni to be the same on the right and on the left hand side of the diagram. Concerning 
the produced pairs, a quick inspection to the kinematics shows that the component Q+ comes 
from and the component Q_ comes from a_, which, neglecting terms of order l/\/s, imply 
ii+ = l'i+ and ai_ = a'i_. In terms of the fractional momenta Xj = x'^^ Zi = z\ and according 
with the definitions 



Using the already defined V'3(Z) V'3(a) and analogously for the factors depending on /'j, a'j, we 
perform a Fourier transform on the transverse momenta; ^ is the conjugated of I and h is the 
conjugated of a. 

The subsequent integrations involve the nuclear variables. The longitudinal variables Z are 
common to both sides of the cut diagram, the transverse variables are different. In the case of 
the Deuteron one has to integrate over one spectator (which is on shell) so we have integrations 
in dZi/Zi and in dNi_i, note that A^i_|_ = — A^2± and Zi = 2 — Z2. In the Tritium i^He) case 
we have two transverse variables and the longitudinal integration, which may be expressed as: 
/ 5{Z\ + Z2 + Z3 — 2>)dZidZ2/ Z2dZ^/ Z^. In conclusion in both cases only one nuclear variable 
survives, and we need to perform the Fourier transform of the transverse components: we call 
the corresponding coordinates B,B'. 

The integrations which take care of the conservation conditions give two kinds of results: the 
diagonalization in the impact parameters given as: B = B' , /3i = f3'-, hi = h[ and the geometrical 
conditions : /3i — /32 = ^1 — ^2, /32 — /Sa = 62 — &3, which imply also /Ss — /3i = 63 — 

More in detail, the derivation of the condition on the transverse variables is obtained as follows 
(to simplify the notation, the transverse index, like l±, is everywhere understood): 

The Fourier transform of the wave function is il^3{li,l2,h) = J ''PsiPi) e'^p[iJ2hPi] TldP and 
analogously for the functions in a and for the conjugated. The conservation relations involving 
the produced pairs, integrated over the final states, give the two-dimensional constraints 



The conservation between the incoming momenta and the momenta of the remnants gives a 
factor S(Y1 ^ ~ ^') and similar for the Oj, the sum is redundant, it is already contained in the 



V'3,m(o) 



1 



(19) 



2iV_2 ZiZ2a3']_ + Z2Z3ai']_ + z-^zxa-;^^ -)- \M\I(\ -Yli^i)- 'm%]ziZ2Z3 ' 



Y[dl+da- = {L+N_fYldxdz . 
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previous relations, the difference gives a new condition S{Y^ / — ^ — ^ a + ^ a') . In exponential 
form: 



(2vr)2 

The integrations over the internal variables I, I' , a, a' give; 

6{(3i + ei-C), 6{/3'i + 9i-C), 6{bi + 9i-C), 6{b'i + 9i-C) 

which in turn implies: /3j = bi = b[. One thus obtains the diagonalization in the impact 

parameter. Moreover one is left with 6{l3i — bi — 2C), which represents a geometrical constraint: 
the difference /3j — bi is independent of the index i. 

In analogy with the previous definitions |55] [58]. the three-body densities are defined reab- 
sorbing partially the factors L-|_, 

1 / ~ X\ 

T{xi,X2,X2,;[ii,(32,fiz) = .a / iV^la^i, 2^2, 2:3; /32, /Ss)^^ — L+'^dM} 

2(27r)y J I - xi - X2 - X'i 

T{zi,Z2,zs;bi,b2,b3) = — J-g / 1^^(^1,22,^3; 61,^2,53)1':^ ^-^^ NJdMf (20) 

2 [Itt)^ J 1- Zl- Z2- Z3 

In this way the densities T are such that they are neither vanishing nor growing indefinitely 
when L+, N- — t- 00. The last expression could be recast in terms of the external variables Xi as 

T{xi/Z,X2/Z,X3/Z;bi,b2,b3) (21) 

= Z2— 1- / \4,ixi/Z,X2/Z,xs/Z-M,b2M)\ \ - D.^dM! 

2»(27r)y J Z - xi - X2 - X3 

The triple scattering cross section on a Deuteron, which describes one of the bound nucleons 
suffering three hard collisions, while the other is a spectator, it is hence given by: 

= (2^ y'r(.i,X2,X3;/3i,/32,/33)r(x;/Z,4/Z,x'3/Z;6i,62,53)|^^^ 

|§d(Z; B)\'^Z-^dB dbi db2 dbs dpi d^ d^z 6{bi - 62 - /3i + p2)S{bi - 63 - /3i + Ps) 
dxidx2dx^dxidx2dx^ dZ dQi dVt2 dVt^ . (22) 

When only a single target nucleon interacts with large momentum exchange, nuclear dynamics 
thus takes completely care of the difference between Deuteron and Tritium (or '^He). There 
are some minor differences: the total four momentum of the interacting nucleons con have a 
transverse component. 

3.2. Two different target nucleons interact with large transverse momentum exchange 

3.2.1. General features and diagonal terms 



As already seen, when two or more target nucleons interact with large transverse momentum 
exchange, nuclear and partonic dynamics are interconnected. The presence of the nuclear wave 
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function induces in fact the presence of two kinds of contributions: a "diagonal term" and 
a number of nondiagonal or "interference" terms. The diagonal discontinuity for the triple 
scattering is shown in fig. 6 and its analytical expression is: 



X ri(/iai ^ qiq'i)T2{l2a2 ^ q2q2)T3{l3a3 ^ 93^3) 
X ri*(/'ia'i ^ qWi)n{l2a2 ^ q2q2)T3* {l' 3a' 3 ^ 9393) 

^d{Ni,N2) ^D*{N'i,N'2) 

^ [Nf - w?] [iV| - TV?] [N'l - m?] [N'l - n?\ 

X 8{L-^l- F4)5{Ni - as - Fi)6{N2 - ai - as - F2) 

X (5(L - ^ /' - F4)6iN\ - a's - Fi)S{N'2 - a\ - a' 2 - F2) 

X 5{h +ai- Qi)S{l'i + a'l - Qi)5{l2 + 02 - Q2)S{1'2 + a' 2 - Q2) 

X <5(/3 + aa - Q3M'3 + a' 3 - Q3)5{D -Ni- N2)S{D - N\ - N'2) 

X Yld{^/8)dada'dldl'dFdNdN'dQi (23) 

with the mass shell condition = M^. 

In the case of a Tritium or ^He (with four-momentum T) the corresponding discontinuities 
Disc H(2'1.o) are obtained from Disc A^'^'^'^ by substituting the factor 

^.(iV.iV.) ^n*{N',^N'2) _s^n-N,-N2)SiD-N'.-N'2) (24) 



[Nf - m2][Ar| - m2] [N'l - m'^][N'l - m?] 



with 



^T(jVl,A^2,iV3) ^T*{N\,N'2,N^) 

[Nl - m?] [Ar| - m2] [A^32 - w?] [N'l - m?] [N'l - m?] [N'l - m?] 

X 5{T -N1-N2- N3)S{T - N'l - N'2 - Ns)S{Ni - m 



(25) 
2^ 



The integration runs over dNi, with i = 1,2, 3, and the mass-shell condition is: = M|.. 

Prom the side of the incoming proton the situation and the subsequent manipulations are the 
same as in the previous case, they give rise to the factor -03. From the side of the nucleus wc 
find different structures: the one-body and two-body parton vertices and the singularities which 
put a nucleon on mass shell or the partons on mass shell. 

Here we need the one-parton and the two-parton wave functions, they were already defined 
and an explicit form was given in eq, (5). The longitudinal integration of the remnants is 
performed as in the previous case. We have now three F± and we get: 

l/Fi_ = l/[N.{l-zi-Z2)] l/F2_ = l/[iV_(l-^3)] l/F4+ = l/[L+(l-^x)] . 

Concerning the nuclear variables: since the binding energy is small, the most important 
singularities are those corresponding to the nucleons' mass shell condition. For the Deuteron one 
must thus evaluate l/{Nf — m?) with the other propagator on mass shell, i.e. N2 = m?, with 
A'^i -|- N2 = D. One obtains the expression of ^ i:){N-)/N- as given in Eq (6). For the Tritium 
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Figure 6: Different contributions to the triple parton scattering in p D interactions. Both target nucleons 
interact with large transverse nionientum exchange 

or ^He we have Ni + N2 + N3 = T and, in addition, the mass-shell condition for the spectator, 
we get the expression in eq. (7). The one- two-three parton densities together with their Fourier 
transform in the transverse plain T{xi; hi) have been already defined and discussed. 

Neglecting terms of order the conservation of the large components gives again 



li+ = CL\_ = a'i_ and the complete expression is brought into a diagonal form by Fourier 
transforming the transverse variables. In this case, the geometrical relations are however differ- 
ent. One finds: /32 — /^s = 62 — ^3 — ^1 + /3i — /32 = &i — &2 + ^1 — -B2. The corresponding 
configuration in transverse space is shown in Fig 7a). 

The diagonal contribution to the triple parton scattering cross sections on Deuteron and on 
Tritium or ^i?e, when one nucleon interacts twice and another once, setting B = B\ — B2 and 
B' = Bi — B3, are thus expressed as: 



pD\ 
^3,2 \d 



r(xi, X2, X3; /3i, /32, (33)T{zi,Z2; bi,b2)T{z:i; 63) 

b2 + (3i 



da da da 



(2^)3 

\^d{Z; B)fdBdPid^2dP3 dbidb2db3 6{bi 
dxidx2dx3 dzidz2dz3 dQid^}2dfl3 dZjZ^ 



dVL\ dQ2 
P2) d{bi - 63 - /3i + /33 - 5) 
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-fiy^ f r(x„x„x3;A.A.A)r(.,...;6„fc)rfe;63)||||Ji (26) 

5{bi — 62 + /3i — P2) 5{bi — 63 — /5i + /33 — 5) dxidx2dxz dzidz2dz^ 
dnidn2dn3 5{Zi + Z2 + Z3- 3)dZidZ2dZ3/{ZfZ2) . 



3.2.2. Interference terms 



As it appears from the graphs in Fig. 6, there are three kind of interference terms which differ 
one another for the different relations between the partons and the parent nucleon. 

This difference takes, therefore, the form of difference in the (5— functions: the discontinuities 
corresponding to three terms are expHcitly given for the case of the deuteron: they can be 
summarized in the following form: 



\^n (2^)35 ] q^q iqiq^q „2„2 „2 ^,2 ^,2^,2 

X Ti(/iai qiq'i)T2{l2a2 9292)^3(^303 '?3'?3) 
X ri*(ria'i qWi)n{ha2 ^ q2q2)n* {I' 3a 3 ^ QsQa) 
^d{Ni,N2) ^D*{N'i,N'2) 

^ [iVf _ ^2] [jv| _ [N'l - m?] [N'l - m?] 

X 5{L-^l- Fi)S{L - ^ r - F4) X r]{Ni, aj) 

X S{h +ai- Qi)d{l'i + a I - Qi)S{l2 + 02 - Q2)5{1'2 + a 2 - Q2) 

X 5(13 + 03 - Q3)Sil'3 + fl'3 - Q3)SiD - iVi - A^2)(^(^ - A^'i - A^'2) 

X Yld^dada'dldl'dFdNdN'dQ (27) 

The three relevant reahzations of the factor r), as it can be seen from the graphs, are: 



r?2,i(^i, = <^(^i - «3 - Fi)S{N2 -ai-a2- F2)5(7V{ - ^ - Fi),5(iV^ - a[ -a'3- F2) 
V2,2im, aj) = 5{Ni -as- Fi)6{N2 - ai - 02 - F2)5{N[ - a[ - 4 - Fi)S{N^ - 4 - F2) 
md^i, aj) = 6{Ni -as- Fi)6{N2 - ai - 02 - F2)d{Ni - a[ - a'2 - Fi)(5(Ar^ - 4 - F2) (28) 

The corresponding expression for the case of the Tritium are obtained by the same substitu- 
tions that were used in the diagonal term. 

The configuration produced by the factors ri2j are similar to the configuration described by 
the crossed diagram in the double scattering, however the factors are more strictly interlocked 
so that it is necesssary to introduce other auxiliary terms 



. 1 X1X2X3 

h,i {xi,X2,X3;Zi,Z2;bi,b2,b3;B) = — — — = — 

4:{2ny (Zi - xi){Z2 - X2 - X3) 

X J tpl,Miixi/Zi;bi)'4)2,M2{x2/Z2,X3/Z2;b2,b3) 

X r2,MA^i/Z'2,X3/Z!2;b3,b2 + B)ri,MA3;2/Z[-,B-bi)dMUMi (29) 
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The nuclear factors have formally the same expression as in the diagonal case, but the values 
of Z are different on the two sides, while keeping the constraints Zi + Z2 = Z'^^ -\- Z2 = 2 . We 
have in fact Zi — xi = Z'^^ — X2, Z^ — = Z'2 — x\. 




Ai=|3i-|32 

A2=P3-P2 



B B 



c) d) 



Figure 7: Configurations in transverse space of the four amplitudes in Fig. 6, in the right hand side of tlie 
cut 



H^2,2 (xi,X2,a;3;Zi,Z2;6i,62,63;5) ^ ^\^2Xz 



4(27r)9 (Zi -xi)(Z2-:e2-X3) 

X y V'i,Afi(Si/-^i;&i)V'2,M2 (^2/^2, ^3/^2; ^2,^*3) 

X xs/Zl; 62 + B, hiWiM^i^^/Z'^, bs)dMf dM^ (30) 

Here beyond the the constraints Zi + Z2 = + = 2 we have the relations Zi = Z'^ — 

Xs, Z2 = Z'2 + X2,. 



W^2,3 (xi,X2,X3;Zi,Z2;6i,62,63;-B) ^i^^x^ 



4(2^)9 (Zi-Xi)(Z2-X2-X3) 

X j V'i,A/i(Si/2'i;6i)-02,Af2 (^2/^2,^3/^2; 63) 

r2M. {X2/Z[,x^/Z[; b2 + BM + BWiM. (^1/^2; B - h)dM^ dM^ (31) 



X 
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Here beyond the the constraints Zi + Z2 = Z[ + Z2 = 2 we have the relations Zi — xi = 

Z[- X2- Xs, Z2 - X2 - X3 = Z'2 - Xi. 

We saw that the nuclear factors have formally the same expression as in the diagonal part, 
but the values of Z are different on the two sides. This feature was already found in the cross 
diagram for double scattering, and so also the qualitative consideration are of the same kind. 
Precisely the longitudinal variables arc on one side Zi and Z2 whereas on the other side they are 
Z'l and Z'2. In the nonrelativistic conditions of the internal motion, in particular « m?, 
which are the actual conditions in the Deuteron and Tritium the typical width of the nuclear 
wave function of the Deuteron, in dimensionless variables, is of the order {Am? — M^)/(M|,), 
while the differences Z — Z' are of the order of the fractional momenta x. 

Now we see that the cross sections for these particular processes can be obtained from the 
expression of the diagonal term (eq.25) by substituting the factors V{zi, Z2\b\,h2)V{zz\hz) with 
the correspnding W2J term, after solving the constraints which give Z'-^^Z'^ in terms of Z\^Z2- 

In comparing the diagonal term with the interference terms, one can repeat the same qual- 
itative considerations made for the double scattering case, i.e. in the partonic amplitudes the 
denominators Z., Z' could be set equal to 1, while their value is relevant in the nuclear functions; 
the transverse variables are of the order of the hard size of the hadron, so they are relatively 
small as compared with the nuclear size. 



3.3. Three different target nucleons interact with large transverse momentum exchange 



3.S.I. General features and diagonal terms 



This kind of process can evidently happen only with a three-body nucleus (at least), here also 
the presence of the nuclear wave function induces the presence of two kinds of contributions: a 
"diagonal term" and a number of nondiagonal or "interference" terms. 

The diagonal discontinuity for the triple scattering is: 



Disc 



_ 1 f 4>p 4>*p (pi 4>2 4>3 (f>l 03* 

d ~ (27r)38 J qqq iqiqiq a\ 4 4 aq^l ai 

X Txilxax qiq[)T2{l2a2 -> 92^2)^3(^303 Qsgs) 

X ri*(ria'i ^ qiq[)nil2aa ^ q2q2W{l'3a's ^ 53^3) 

^ <^D{Ni,N2,N-i) ^D*{N\,N'2),N'-s) 

^ [Nf - m2] [iV| - m2] [Nl - m^] [N'l - m"^] [N'l - m^] [Nq - n?] 

X 5{L-Y,l- Fi)5{Ni -ai- Fi)5{N2 - 03 - F2)5{N3 - 03 - F3) 

X 5(L - ^ r - Fi)5{N'i - all - Fi)5{N'2 - a 2 - F2)b{N' j, - a'3 - F3) 

X 5(^1 + ai - Qi)b{l\ + a'l - Qx)b{l2 + 02 - Q2)K^ 2 + a'2 - Q2) 

X S{h + aa - Q3)(5(r3 + a'3 - Q3)<5(r - TVi - ATs - Nz)b{T - A^'i - N'2 - N'3) 

X Yldi^/8)dada'dldl'dFdNdN'dQ (32) 

and the mass shell condition is = M^. 

From the projectile side we have the three-parton densities, which have been already defined 
and used. On the Tritium side we must use the full three-body structure of the nuclear wave 
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c) 



Figure 8: Different contributions to triple parton scattering in p T or p ^He interactions. All target 
nucleons interact with large transverse momentum exchange 



function: the conservation iVj 
performed may be written as: 



a-i + Fi gives a} = {Ni - Fi] 



Mi. 



ai\ then the integral to be 



ya= I dN^+dN2+dNs+s{ iv,+ - r+) Yl 

1=1,2,3 j=l,2,3 



As in the previous cases there are two kinds of singularities in the integrand, one kind puts a 
nucleon on mass shell, the other one puts on mass shell the parton, and here also the singularities 
putting the nucleons on mass shell are the most important ones. So we approximate as: 



ya 



1 



h-h-h 



n 



.=1,2,3 K< - N-F+\a,_ - (iV-ai),] 



mx'/A^-, (33) 



When considering the complete cut graph we find that the relation between l,a,Q yields 
again the equalities /+ = a_ = a'_, whereas the equality of the remnants yields = A^'_, 
these can be converted into the fractional momenta Z, with the constraint Zi + Z2 + Z-^ = 3. 
Here also it is convenient to go from transverse momenta to transverse coordinates. 

The factors ip and T are defined as before. Since in ya the dependences on a± and on N± 
are interlocked, the Fourier transform must be performed with respect to both sets of transverse 
variables a± and A^_l. 
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Figure 9: Configurations in transverse space of the three amplitudes in Fig.8, in the right hand side of 
the cut 



As before /3j are conjugated to li, hi are conjugated to and Bi are conjugated to iVj. The 
conservations, that will be again conveniently expressed in exponential form, are: li + m = 
I'i + a'i, = "^I'i, Ni + Qi = N'l + a'i, — = ^- 'T^^ integrations over the 

transverse momenta yield the equalities Pi = f3'-, bi = b[, Bi = B'- and the geometrical conditions 
bi + Bi — /3i = const that can be expressed also as: bi + Bi — f3i = bj + Bj — (3j; i ^ j 

We can write, for the diagonal term: 



^asL = f ^{xi,X2,X3; I3i, (32, P3)'r{zi;bi)T{z2;b2)T{z3;b3) — — — 



\^TiZi;B, B)\'^dBdBdf3idP2dp3 dbidb2db-i 
5{bi - 63 - /3i + /33 + 5{bi -b2- f3i + 1^2 + B) dxidx2dx3 dzidz2dz3 
dnidn2dn3 5{Zi + Z2 + Z3- 3)dZidZ2dZ3/{ZiZ2f . (34) 

In the vertex function only the differences Bi — Bj are relevant, therefore the integration 
variables B,B represent a pair of these differences, e.g. B = Bi — B2, B^ = Bi — B^^. 
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3.3.2. Interference terms 



The discontinuity for the interference terms in triple scattering can be written, as in the 
previous chapter: 



n- 1-1)1 = _*_ f ^* 

(2;r)38 J qqq iqiqiq aj 4 a'la'l 4 

X ri(/iai qiq[)T2{l2a2 9292)^3(^303 qsq's) 

X Ti*{l'iai qiqi)T2{l2a2 924)^3* (^'3«'3 ^ ^s^s) 

^d{Ni,N2,N^) ^D*{N'i,N'2),N'^) 

[Nl - m?] [iV| - m2] m - m2] [N'l - \N<1 - tv?] [N'l - m?] 
X 5{L-Y,l- F^)5{L -^r-F^)x ^{N,, 
X 6{li +ai- Qi)5{l'i + a I - Qi)5{l2 + 02 - Q2)^{1'2 + a 2 - Q2) 
X 6ih + «3 - QsWs + a'3 - QsHT -N1-N2- N3)6{T - N\ - N'2 - iV's) 
X W^d{9./%)dada'dldl'dFdNdN'dQi (35) 

Here we find two essentially different realizations of the factors r/ 





= 5{Ni 


- Ol 


- Fi)6{N2 


- a2 


- F2)6{N3 


- as 


-F3) 




5{N[ 


-«3 


-Fi)6{N^ 


- a'2 


- F2)6{N!, 




-F3) 




= S{Ni 


- fll 


- Fi)5{N2 


- a2 


- F2)6{N3 


- as 


-F3) 




6{N[ 


/ 

- «2 


-Fi)6{N^ 


- as 


- F2)6{N!, 


— ai 


-F3) 



and this leads to the definition of two more auxiliary functions W. 



WsA {xi,X2, X3; Zi,Z2, Z3; bi, 62, 63; B, B) 



2:1X2X3 



4(27r)9(Zi-Xi)(Z2-X2-X3) 
X V'l,Mi(^l/-^i;&l)V'l,M2(^2/^2; &2)V'l,M3(:^3/-^3;^'3) 

X V^t,M3(^l/^i;&l+^)^t,M2(S2/^2;?'2)V^t,Mi(^3/^3;^3-S) (37) 

Here beyond the the constraints Zi + Z2 + = Z'^ + Z2 + Z'^ = 3 we have the relations 
Zi-Xi = Z[- X3, Z2 = Z3 - X3 = Z3 - Xi. 



W^3,2 (xi,X2,X3;Zi,Z2,Z3;6i,62,^s;5,^) = ttttvjT^^ w f ^ 

4(27r)y (Zi - xi)(Z2 - X2 - X3) 

X -01 ,Mi (x 1 /Zi ; 61 ) ^1 (^2 7-^2 ; &2 ) -01 ,M3 (^3 /^3 ; ^3 ) 

X i^lu^ixi/Zi-M - Ai)V'I,M2 (2^2/^2; &2 - ^2)i'lKh{x^/Z-yM + A2 + Ai) (38) 
Here beyond the the constraints Zi + Z2 + Zj, = Z'^ + Z2 + Z'.^ = 3 we have the relations 

Zi - Xi = Z[ - X2, Z2 - X2 = ^2 ~ ^3 -^3 - ^3 = -^3 - ^1- 

Finally the cross sections for these particular processes can be obtained from the expression 



of the diagonal term, Eq.(77), by substituting the factors T{zi]bi)T{z2;b2)T{z^]b^) with the 
correspnding W^j term, after solving the constraints which give Z3 in terms of ^1,^2,^3. 

The transverse configurations are shown in Fig. 9. 
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Clearly the qualitative considerations made previously for ratio between the diagonal and the 
interference term hold also here, since they depend on the existence of two scales (hadronic and 
nuclear) playing always the same role. 



4. SIMPLEST ESTIMATES OF THE DOMINANT CONTRIBUTIONS 



The non-perturbative component of MPI in pA collisions is characterized by the hadronic 
and by the nuclear scale. In MPI, the relevant hadronic scale is the transverse dimension R of 
the generalized parton distributions, which is smaller as compared with the hadron radius and 
it may be roughly a factor four smaller as compared with the radii of D, and ^He. Even 
with light nuclei, one may thus obtain a simplest estimate of the dominant contributions to the 
MPI cross sections by neglecting the hadronic scale when compared to the nuclear scale. In the 
same spirit one may obtain a further simplification by evaluating the integrals on the fractional 
momenta of the bound nucleons, Zj, by keeping into account the dependence on Zi only in the 
nuclear wave function and replacing Zj with 1 anywhere else. 



4.1. Double scattering 

The double parton scattering cross section, for p D and of p (or p ^He) collisions, are 
given by the sum of two contributions, where only a single bound nucleon and where two bound 
nucleons participate in the hard interaction: 

(T2 (^Xi,Xi) = ]^ (3^1, 3; j) + (72 2 (s^i ) 3J j ) (^9) 

With the simplifying assumptions above, the contributions where only a single bound nucleon 
participates are given by 

(7^J ~ 2 (7D and a^'^^ ~ 3 (7£, (40) 

where (7j) is the double parton scattering inclusive cross section on a isolated nucleon. One can 
make the positions: 



r(xi,X2; A,/32) = K^^^^^G{xi)G{x2)fxi,x2{l^i,h) 

with J U,,xM,P2)dPidP2 = l (41) 

where G{x) is the one body inclusive parton distribution, such that G{x) = jT{x;b)db. In the 
case of identical interactions, one may thus express the double parton scattering cross section on 
a isolated nucleon as 



CrDixl,Xi,X2,X2) = ]^Kx^,x2 Kxi ,X2 J fxi ,X2 {Pi > 1^2 ) fxi ,S 2 (^1 > ^2 (a^l , ) (75 {X2 , X2 ) 
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X 5{^i - /32 - foi + b2)d(3idP2dbidb2 (42) 

where erg = J G{x)cr{x, x')G{x')dxdx' is the single parton scattering inclusive cross section on a 
isolated nucleon. The effective cross section, namely the accessible experimental information in 
nucleon-nucleon collisions, is thus given by 



\ — = K^,,^,Ks,^s, [ /xi,x2(/3i,/3i - Ai)/^,,^,(6i,6i - Ai)dl3idbidAi (43) 

aeff{xi,Xl,X2,X2) J 

where Ai = /3i — (32- In p D collisions, when two nucleons participate to the hard interaction, 
one has contributions from a diagonal and from a off diagonal term. The dominant contribution 
to the diagonal term is 

'^2^L(^i'^i'^2,^2) = jT{xi,X2;f3i,P2)^ixi,xi)a{x2,X2)T{xi/Z;bi)T{x2/i2-Z);b2) 

X \^D{Z;B)\^dZ/Z^ dB dbidb2 dpi dP2S{B -bi + b2- Ai) 

~ K^^^^^as{xi,xi)as{x2,X2)'lDiO) (44) 

where the general form of Xz)(x) is: 



Id{x) ^ / ^D{Z;0)n{Z';0)^^6{Z - Z' -x) (45) 

The range of the variables bi and Pj, defined by the partonic distributions F, is much narrower 



as compared with the nuclear range of B. The (^-function in Eq. (44) thus forces, in ^/j, ~ 0. 

Notice that, differently from the case of double parton scattering in nucleon-nucleon collisions, 
the scale factor in (^2^1^ is given by the value of I/)(0), which is determined by the radius of the 
Deuteron, while the cross section is proportional to Kxj^^x2i which gives the partonic correlation 
in fractional momenta. As already noticed in [58], crf?]^ depends thus weakly on the correlation 
between partons in the transverse coordinates and, on the contrary, it may provide a rather 
direct information on the size of Kxj^^x2- 

The contribution of the interference term to the cross section is 

cr^^l (xi,xi,X2,X2) = -i-^-TK I T{xi,X2;Pi, P2)^{xi,xi)a{x2,X2)Wi{Z,Z';xi,X2;bi,b2,B) 



X ^D{Z;B)^}y{Z';B)[ZZ'] ^dB dbidb2dpidp2dZdZ' 



(27r)3 

X 6{B -bi + b2- Ai)6{Z - Z' -X1 + X2) (46) 
By neglecting the hadron scale with respect to the nuclear scale and keeping Z ^ 1 only in 



the Deuteron wave function, the integrations in bi and 62 in Eq.(46) are 



J Wi{l,l;xi,X2;bi,b2,bi - 62 + Ai)d6i db2 



dMldM^dbi db2- — — V'Afi(xi; biWM^{x2; 62 - Ai) 



4(27r)6 7 ^ ^ -(i_5;,)(i_s2) 
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xtj;M2{xi;b2)ip*Mj^{x2;bi + Ai) 
= ^(xi,X2;Ai)^(a;2,a;i;-Ai) (47) 

where the generahzed parton distributions H have been introduced: 

#(xi,X2;Ai) = J Hil,l;xi,X2;bi,b2)S{h-b2- Ai)dh db2 (48) 
Notice that the normahzation is H{x,x;0) = G{x). One may thus define 

C,{x„x„X2,X2) = / f..M^i)H{xuX2;A^^^^^^ dA, 

G{xi)G{x2) 

which is dimensionless and weakly dependent of xi, X2 as compared to Id(xi — X2), since 
Ci originates from the partonic structure of the hadron, while X^) originates from the nuclear 
structure. The contribution of the interference term to the cross section may thus be expressed 
as 

(^2^\Axi,Xl,X2,X2) ^ Kx^^x2<^s{xi,Xi)as{x2,X2) Ci{xi, Xi, X2, X2)Id{xi - X2) (50) 



Notice that both (^2^1^ and cr2^\- depend linearly on Kx-^,X2 and both terms are proportional 
to the inverse of the square of the Deuteron radius, the latter term through the nuclear off 
diagonal factor X£){xi — X2), which induces a much stronger dependence of o'2^\. on xi — X2 as 
compared with crf^]^. 

In the case of double parton interactions in p or p ^He collisions, with two target nu- 
cleons taking part to the hard interaction, one obtains the same expressions, for the dominant 
contributions to the cross sections, as in the case of p D interactions. The only difference is 
in the multiplicity factors and in the terms loiO) and l£){xi — X2), which are replaced by the 
corresponding quantities with or ^He, actually IxiO) and Ixixi — ^2). 

The leading contributions to the double parton scattering cross sections in p D and p ^H, 
p ^He are thus given by 



^2 ("^j' "^i) — ^2 1 ('^*' "^i) ~^ ^2 2 l<i('^i' •^i) ~^ ^2 2 I ^i) 

~ 2 aD{xi,Xi) + Kxj^,x2C^sixi,xi)asix2,X2)[lDiO) + Ci{xi, Xi)lDixi - X2)] 

) (51) 
~ 3 aDixi,Xi) + 3 Kxi^x2<^s{xi,xi)asix2,X2)[lT{0) + Ci{xi,Xi)lTixi - X2)] 

The contributions and af^ are well approximated by 2 ctd and by 3 ao- The actual 
values can be evaluated with great accuracy, once the double parton scattering cross sections in 
pp and in pn collisions are known as a function of fractional momenta. Also the off diagonal 
nuclear terms X/)(0), Id{xi — X2), 2r(0) and Ixixi — X2) can be evaluated very accurately. By 
measuring the double parton scattering cross sections in p D and p (or p '^He) one may thus 
obtain accurate estimates of the differences (72^ — o"!^ and o"!^ — and, as a consequence, of 
the rations TZd and TZt, defined as 
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nD{Xi,Xi) = T — ~ Aa;i,a;2 lz)(0) + Ci(Xi,Xj)Zz)(Xi - X2) 

crs{xi,xi)a3[X2,X2) 

pT I , - \ _ pT ( - \ 
TZT{Xi,Xi)= — / ~ 3 Kxi,x2\^TiO) + Cl{xi,Xi)lT{xi - X2)\ (52) 

(7s{xi,xi)as{x2,X2) 

One should point out that the contribution of the interference term is not always present. 
As an example, in the case of production of a l^/'+jets, through double parton collisions, the 



interference term is absent. In such a case, the second term in square brackets in (52) is missing 



and Eq. (52) allows a direct estimate of Kx^^x2: namely of the importance of the correlations in 
X in the double parton distributions. When the interference term is present, the cross section 
depends on the additional unknown quantity, Ci{xi,Xi), which multiplies the nuclear overlap 
integrals Xd{xi — X2) or It{xi — X2), which originate the main dependence of the cross section 
on X = xi — X2. By measuring the cross section at different values of X one may construct the 
fraction 

^^^^^ - K.,„C.i^....)^^ (53) 
By studying the dependence of TZd and of TZ^ on X = xi — X2, with the help of Eq.s (|52|) 



and (53), one may obtain information both on Ci{xi,Xi) and on Kxj^,x2- Notice that the values 



of AX needed in Eq.(53) may not be too small. The natural scale of X = xi — X2 is in fact 
{Eb /TnY^'^ ~ 5 X 10^^, where Eb is the nuclear binding energy. 

The indication on the value of Kx^,x2^ together with the measure of the effective 
cross section in nucleon-nucleon collisions, allow to obtain an indication on the value of 



the integral / fx^,x2{^i)fx^,x2{^i)d^i (cfr Eq. ^). With the help of Eq.s ([52]) and 
( |53p , one may estimate Ci(xj,Xj) and obtain an indication on the value of the integral 
/ fxx,x2{^i)^{^i-i^2]A.i)H{x2TXi;—/S.i) dAi (cfr Eq. ([49^). The information on the two inte- 
grals will provide important constraints on the correlation length between partons in transverse 
space, which is explicit in fx^,x2{^i)- By measuring the double parton scattering cross section 
Yo. p D and p (or p ^He) one may thus learn both on correlations between partons in frac- 
tional momenta, through the factor Kxi^x2^ ^iid 011 correlations between partons in the transverse 
coordinates. 



4.2. Triple scattering 

As in the case of double parton scattering, the triple parton scattering cross section \n p D 
and of p or p ^He collisions can be written as 

C3 (^Xi,Xi) = (T3 -|^(xj, Xj) -|- £73 2(2^1) -|- (T3 3(xj, Xj) (54) 

With the simplifying assumptions discussed in the previous chapter, the dominant contribu- 
tion to the terms where only a single bound nucleon undergoes a triple parton interaction is 
given by 
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f^3? - 2 (7T, af J ~ 3 (TT (55) 

where ut is the triple parton scattering inclusive cross section on a isolated nucleon. By making 
the positions: 



{Pl,P2,P3) (56) 

with 



/ 



fxi ,X2 ,X3 {Pi,(32,P3)dPidl32dl33 = l (57) 



one may express, in the case of identical interactions, the triple parton scattering cross section 
on a isolated nucleon as 

1 f 

(^T{xi,Xi,X2,X2,X'i,X'i) = — -K3;i,a;2,X3-^Si,S2,23 / fxi,X2,X3 ,X2,X3 (bi) (58) 

X 5(^1 -P2-h + b2)S{Pi -P3-h + b3)Ylas{xi,Xi)dPidbi 

4-2.1. Two different target nucleons interact with large transverse momentum exchange 

The contribution to the triple parton scattering cross section, where two target nucleons 
undergo hard interactions, is a process of 0(l/(5^i?^)), to be compared with triple scattering on 
a single nuceon, which is of ©(l/i?*^), R and S are the hadronic and nuclear scales. In the case of 
p D collisions, the different contributions to the cross section are summarized by the expression: 



= (2^-^^7^(^^'^^'^^'^^'^^'^3)d|||d| 

X W2,j{xx,X2, xs; Zi, Z2; bi, b2, 63; B)^d{Zi, Z2; B)^}j{Z[, Z'^, B) 
X dB dZidZ[dZ2dZ'2/{ZiZ[)Y[dxidxidbidPidnidZi 

X S{Zi + Z2- 2)SiZ[ + Z'2- 2)5{Zi -Z[- Xj)5{Z2 - Z'^ + Xj) 

X S{b2 -bi + Ai)S{B - 61 + 63 - A2) (59) 

where Ai = /3i — /?2 and A2 = /?3 — while the index j corresponds to the diagonal case, 
when j = 0, and to the the three different interference terms, when j = 1,2,3. By neglecting 
the hadronic scale R as compared to the nuclear scale S, in the diagonal case one obtains 

W2,oixi,X2,X3; Zi, Z2;bi,b2,b3; B) = ^T{zi, Z2;bi,b2)T{zs;b3) (60) 
The quantity Xj assumes the following values: 
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Xo = 0, Xi=Xi-X2, X2 = -X2, X2=Xi-X2-X3 (61) 

and the multiplicity factors Afj are: A/o = 2, M = 4, J\f2 = 4, A/s = 2. 

Analogously to the case previously discussed, the dominant contributions may be estimated 

by 



a 



pD 
3,2 



1^. ~ Afj J T{xi,X2,X3;(3i,l32,/33)a{xi,xi)&{x2,X2)a{x3,X3) 

X W2j{xi,X2,X3;l,l;bi,b2,b3;bi - 63 + A2)X£i(Xj) 

X (5(62 - &i + Ai) JJ dxidxidbi dPi (62) 

i 

where all effects of the Deuteron wave function are summarized in the terms Xu^Xj): 

^D{Xi) = J^,j ^D{Z;0)n{Z';0)S{Z - Z' - X,)^^ (63) 
already defined in Eq.45. By introducing 



/ fx^,x2,x3{Pi)W2,j{xi; l, l;bi;bi-b3 + A2 )5{b2 - 61 + Ai) U- dhP, 

G(xi)G(x2)G(x3 

one obtains 



(^s^lji^h^i) - MjKxias{xi,xi)as{x2,X2)as{x3,X3)C2,j{xi,Xi)lD{Xj) (65) 
The expression of the diagonal contribution is 

<^3^\o{Xi>^i) - ^K:cuX2,X3Ksi^X2O^s(.Xl,Xl)as{x2,X2)as{x3,X3)gx^s^^ID{0) (66) 

where 

GxiXk = J fxi,x2,x3{Pi,P2,/33)fxix2{bi,b2)S{b2-bi + Ai)Y[d/3idbk (67) 
As in the case of double collisions one may introduce the ratio 

and, for j > 0, construct the fraction 

. KMi.,,.,)^ (69) 
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Differently witli respect to tlie case of the double collisions, Eq.s(68) and (69) do not allow 
to disentangle K^^ from C2j{xi, Xi). Disentangling the effects of longitudinal and transverse 
correlations was possible in the case of double collisions because, in that case, the dominant 
contribution to the diagonal term depends only on Kxi^x2- Iii the actual case, on the contrary, 
C2fl{xi, Xi) is proportional to the product Kx^QxiXk- By studying triple scattering on Deuteron 
one may only obtain an estimate of the products KxiC2j{xi,Xi). To gain further insight into 
longitudinal and transverse three-body correlations one needs additional information, which can 
be provided by triple parton interactions in collisions of protons with or with ^He. 

When two target nucleons participate to the hard interaction, after integrating on the spec- 
tator nucleon, p "^H (or p ^He) give very similar results as p D collisions. The two dominant 
contributions to the diagonal term differ in fact only by an overall multiplicity factor (which is 
actually 3) and in the factors Xj'(Xj), which substitute the factors lui^j)- Analogously to the 
case of L), IxiXj) represent the square of the (or ^He) wave function in the mixed represen- 
tation, integrated in the fractional momenta Zi with the constraints given in Eq. (59), and in the 
relative transverse distance B' , while the transverse distance B has been set equal to zero. One 
thus obtains the relation 



^3,2 



Eq. ((to]) is a consequence of Eq.(|68|, which holds in the limit R^jS^, 0. If St is the 
(or ^He) radius, the two terms in Eq. (70) are thus of 0^1/(i2^S'|i)^ and the relation is exact 
in the limit / S'j^ — t- 0. Finite values of R contribute, in the left hand side, with terms of 
0{l/{R^S'^) X R^/Sl^ and with terms of o(l/{R^Sl) x R^ / S^^ in the right hand side of the 

equation. One may thus estimate that Eq. ( |70| ) is valid up to terms of O ('l/ S'^ x (1/ S"^ — 1/ S"^) ~ 

{S'jj — S^)/Sj, 



4-2.2. Three different target nucleons interact with large transverse momentum exchange 

In the case of the contribution to the triple parton scattering cross section, cr^'^, where three 
different target nucleons interact with large transverse momentum exchange, in p or p ^He 
collisions one has three different terms, one diagonal and two off-diagonal, which are labeled with 
the index j in the expression here below. As in the previous section, the label j = corresponds 
to the diagonal case. 



■SI. = (^-^'/r^^^'^^'^^'^i'^^'^^)!^^^ 

X W3,j{xi,X2,X3; ZuZ2, Zs; h, 62, h; B, B')^TiZi, Z2, Z-y, B, B') 

X ^*t{Zi,' Z'2,Z'y,B,B')dB dZidZ[dZ2dZ2dZ3dZy{ZiZ[Z2Z2) 

X 6{Zi + Z2 + Z3- 3)6{Z[ + Z'2 + Z'3- 3)5(Zi - Z[ - Yj^i)6{Z2 - + ^^-,2) 

X SiZs -Z's + y,- 3) SiB' + 61 - 63 + Ai) 6ib3 - 62 - A2 - 

X YldxidxidbidPidnidBdB' (71) 

i 
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With the approximations previously discussed one obtains 
2 



=^'=^1^- - (27r) 



3 A/} / r{xi,X2,X3;l3i,P2,l33)&{xi,Xi)a{x2,X2)a{x3,X3) 



X W3j{xi,X2,X3;l,l,l;bi,b2,b3;b3 - b2 - A2, -61 + 63 - Ai) 

X J{Yj,i=i,3) n dxidxidbi dPi (72) 



where 



J{Yj,i=i,3) = J ^t{Zi,Z2,Z3;0,0)^*t{Z[,Z!„Z!,;0,0) 



X dZidZ[dZ2dZ'2dZ3dZ'2 [ZiZjZaZ^ZsZ^] 
X S{Zi + Z2 + Z3- 3)S{Z[ + Z'2 + Z'^-2) 
X 5{Zr-Z[-Yj,r)5{Z2-Z'2 + Y^,2) 



-1/2 



and 



^j=0,i=i,3 : 


^0,1 


= ; 




Yo,2 


= ; 


Yo,3 


= 




l^=i,j=i,3 : 


^1,1 


= Xi 


- X3 




= ; 


^1,3 


= X3 


- Xl 


5^=2,1=1,3 : 


^2,1 


= Xl 


- X2 


l2,2 


= X2- X3 


Y2,3 


= X3 


- Xl 



while the multiplicity factors are Mq = 1, M[ = 3, M[ = 2. Introducing 



C3,j{Xi,Xi) = 

J fxi,X2,X3 

{Pi)W3j{xi;l,l,l;bi;b3 



b2 - A2, -61 + 63 - Ai) H- dbi dPi 



G{xi)G{x2)Gix3) 



one obtains 



C3,oiXi,Xi) = 1 



and one obtains 



(73) 



(74) 



a^'^^{xi,Xi)\. ~ M-Kx^,x2,x3(^s{xi,xi)as{x2,X2)as{x3,X3)C3j{xi,Xi)J{Yj^i=i^3) (75) 
In the case j = one has 



(76) 



•^S.aU^*'^*) - Q^xi,X2,X3'^sixi,Xi)asix2,X2)crsix3,X3)J{0) 

For j = 1 one has 

y W3,i(xj; 1, 1, 1; bi- 63 - 62 - A2, -61 + 63 - Ai) JJ dh 



(77) 
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= Hixi,X3; Ai)G(x2)^(x3, xy, -Ai) (78) 

in such a way that 

^ / - ^ _ / fx,,x2,xs{^i)H{xi,X3; Ai)ff (x3, xi; -Ai)dAi 

G(xi)G(x3) ^ ^ 

where 

(Ai) = y /,,,,„,3(A)<5(Ai -/3i +/32)nrfA (80) 

For j = 2 one has 

y" W3,2(xi; 1, 1, 1; bi; 63 - ^2 - A2, -61 + 63 - Ai) ]J d6i 

= Hixi,X3; Ai + A2)i?(x2, xi; -A2)^(x3, X2; -Ai) (81) 



and 



C3,2(a;i,2;j) = 

/ /xi,a:2,X3(Ai, A2)g(xi,X3; Ai + A2)-H'(x2 , Xj ; -A2)g(x3,X2; -Ai)(jAi(jA2 

G(xi)G(x2)G(x3) 



(82) 



where 



/,,,,,,,3(Ai,A2) = I fx„x2,xMn^i - f^i + f^2)S{A2 + Pi - p3)lldl3, (83) 

As in the case of double parton colhsions, discussed previously, the contribution to the triple 
parton scattering cross section cig"^ is given with good approximation by 3 ax, where cjt is the 
triple parton scattering cross section on a isolated nucleon. Once the triple parton scattering cross 
sections in pp and in pn collisions are known as a function of fractional momenta, the smearing 
effects of the nuclear wave function can be taken into account and o"!"^ can be evaluated with great 
accuracy. Also the off diagonal nuclear terms J7(l^-, 1=1,3) can be evaluated with great accuracy. 



By measuring cjg"^, one may thus obtain an accurate estimate of the difference Ug^ — CTg^. Eq. (70) 
allows to estimate CTg"^. One may thus define 

TlTixi,Xi) = ^ — (84) 

where 



^3 (.•^ii'^i) — C3^]^ (Xj, Xj) + (Tg 2 (2^1 ) ^^j) ~l~ Cr3,3 (^j 1 ^«) 
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'^3,2 (-^i ' -^i) ~ '^3.2(-^^' — '^3,2 ( 



(85) 



and one may thus relate the 'known' quantity TZ'rp to the unknown properties of the hadron 
structure, represented by K^^ and C3j{xi, Xi) : 



T^rpi^Xi, Xi 



Tlj "^3,3 -^i) 

Y[as{xi,Xi) 



K, 



j 



M;C3jiXi,Xi)J{Yj^i=l,3) 



(86) 



Analogously to the case of double parton scattering, an indication on triple correlations 
in fractional momenta and in the transverse coordinates can then be obtained looking at the 
variation of TZ'rp as a function of Yj^i and using the property that Cs^o = 1- 

Notice that Eq.(86 ) is a consequence of Eq. (70), which holds up to terms of O (^{S'j^ — Sp) / . 

Since the right hand side of Eq. (86) is of o(l/Sp^ , one may estimate that the relative correction 
to the dominant terms in Eq. (86) is only of Oi{S'j-) — Sp)/Sp ~ 1/5 )■ Eq. (86) can therefore 



provide only semi-quantitative indication on the size of triple correlations, while a better deter- 
mination requires a dedicated study. 



5. FINAL DISCUSSION 

MPI in pA collisions allow obtaining information on multi-parton correlations, which cannot 
be provided by studying MPI in pp collisions [54]. Relevant features of the simplest case, double 
parton interactions in pD collisions, were pointed out in |58j . In the present paper, we have 
extended the study of MPI in pA collisions to the cases of double and triple parton interactions 
in collisions of protons with D, and '^He, including in the discussion also the effects of 
interference terms. 

Double parton interactions in collisions of protons with D, '^H and '^He are discussed in Sec. 2. 
When only a single nucleon takes part to the hard process (Sub-Sec.2.1), the integrations on the 
relative transverse coordinates of the spectator nucleons are decoupled from all other transverse 
variables and the cross section is the same as measured in nucleon-nucleon collisions; apart from 
the proper multiplicity factor and the smearing corrections in the longitudinal variables, which, 
as discussed in Section 4, are however rather small. The explicit expression of the cross sections 
for pD, p and p ^He collisions are given in Eq.lO. Similar considerations hold in the case of 
triple parton collisions on a single nucleon, which is discussed in Sub-Sec. 3.1. The corresponding 
contribution to the triple parton scattering cross section is given in Eq. 21. Notice that the 
spectator nucleons are on mass shell. As already discussed in |58] . in spite of that, one may still 
claim that final state interactions of the spectators are taken approximately into account. The 
statement is supported by unitarity: If a nucleon is produced on mass shell and undergoes a final 
state interaction with the remnants of another nucleon, final state interaction does not modify 
the inclusive cross section, since the spectators are not observed. If a nucleon is produced off 
mass shell, its virtuality is anyhow rather small and it may not be unreasonable to extend the 
unitarity relation SS^ = 1 to the actual kinematical domain. Unitarity allows hence replacing 
all final state interaction with cut nucleon lines, i.e. with on mass shell nucleons. 
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In Sub Section 2.2, we discuss the case of double parton collisions, with two nucleons taking 
part to the hard process. In addition to the diagonal contribution in Fig. 2a discussed in |58) . 
which leads to the geometrical picture of the interaction in transverse space shown in Fig. 3, one 
has a contribution from the non-diagonal contribution in Fig. 2b. The geometrical picture in 
transverse space, of the corresponding interfering configurations a) and b), is shown in Fig. 4. 
Notice that, in both interfering configurations, the hard interactions are localized in the same 
points and are well separated in transverse space. As a consequence, the argument for the 
suppression of the interference terms in MPI, discussed in |62], does not apply in this case. Dif- 
ferently from pp collisions, the additional degrees of freedom provided by the nucleus, namely 
the possibility of having different nucleons involved in the hard process, can in fact produce the 
same partonic initial state in different ways, which can thus interfere in the process. As discussed 
in sub section 2.2, the contribution of the interference term is important in the region where the 
fractional momenta of the interacting partons are of order y^E's/m, where Eb and m are the nu- 
clear binding energy and nucleon mass respectively. Differently from the diagonal contributions, 
which are dominated by the most probable nuclear configuration, where all nucleon's fractional 
momenta are equal, in the off-diagonal term of Fig. 2b the fractional momenta Zi, Z[ and Z2, 
Z2, of the two nucleons taking part to the hard process, are forced by kinematics to be different: 
Zi — Z[ = —{Z2 — Z2) = xi — X2- Here xi, X2 are the fractional momenta of the two target 
partons undergoing the double collisions. 

As discussed in sub-section 4.1, the dominant contribution to the interference term can be 
expressed in terms of off diagonal parton distributions, Eq.s(47),(48). The contribution of the 
interference term can be singled out by looking at the dependence of the double parton scattering 
cross section on the difference xi — X2 (cfr. Eq.(|53|)). Keeping into account that the scale, which 
characterizes the dependence of the nuclear wave function on Z, is ^J^Ej^frn^ one may roughly 
estimate that, to single out the contribution of the interference term one needs to measure the 
double parton scattering cross section in an interval x\ — X2 ~ 5 x 10-2 with an accuracy 
grater than 10%. The whole discussion assumes that each couple of scattered partons, and the 
resulting observed particles, can be identified as a definite pair, which requires that each couple 
is sufficiently separated in phase space from the other couples. The quantitative amount of this 
separation depends on the detailed properties of the final state. 

As mentioned in sub-section 4.1, by studying the ratios in Eq.(52) and in Eq.(53), using the 
information on double parton interactions in pp collisions and keeping into account that the 
dominant contribution to the diagonal term depends only on K^^ 1 o^^^ may obtain information 
on parton correlations in fractional momenta and, through the overlap integrals in the transverse 
parton coordinates, which characterize (Jeff and the interference term, also on parton correlations 
in the transverse coordinates. When the interference term is absent, as for production of W+ 
jets, the task of estimating parton correlations may be simpler. In that discussed in|58|. 

all information concerning correlations may be obtained directly from Eq. ( 52 ) . 

The case of triple parton interactions has many features similar to the case of double parton 
interactions. The main difference is in the sizably larger number of contributing terms. As al- 
ready pointed out, the contribution where only a single nucleon participates to the hard process 
is well approximated by the cross section on a isolated nucleon, multiplied by the multiplicity of 
target nucleons, while nuclear smearing effects can give only minor corrections. Things become 
complex when two or three nucleons participate to the hard process. The general features of 
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the contributions where two nucleons participate to the hard process is discussed in Sub Section 
3.2 while, in Sub-sub Section 4.2.1, one may find a simphfied estimate of the different terms. 
In the case of two participating nucleons, one finds a diagonal and three different off diagonal 
contributions. Differently from the case of double parton collisions, in the case of triple parton 
collisions on two different target nucleons, the dominant contribution to the diagonal term de- 
pends both on the correlations in fractional momenta, through i^xi,a;2,x3) and on the correlations 
in the transverse coordinates, through the overlap function Gxi,Xk (defined by Eq. (67)). A con- 
sequence is that one cannot disentangle the effects of longitudinal and transverse correlations in 
triple parton collisions by studying the ratios in Eq. ( 68 ) and in Eq. ( 69 ) in p D interactions only. 

The information on triple parton collisions in p D interactions can be, nevertheless, utilized to 
estimate the contribution to triple parton collisions with two participating nucleons, in the case 
of p or p ^He collisions. By measuring the triple parton scattering cross section in p or 
p ^He collisions, one may thus estimate, using Eq. (70), also the contribution to the cross section 
where three nucleons are involved in the hard process and thus figure out the value of the ratio 
TZ'rp in Eq.(86). A relevant feature of cf^, the component of the triple parton scattering cross 
section with three participating target nucleons, is that the leading contribution to the diagonal 
term, as given by Eq. (77), is proportional to i^a;i,x2,a;3 and does not depend on the correlations 
in the transverse coordinates. By studying the dependence of TZ'rp on Yij one can thus obtain 
an estimate both of i^xi,x2,x3 and of the different overlap integrals in the transverse coordinates, 
which characterize the different interference terms. As discussed in the last part of sub-sub 
section 4.2.2, the uncertainties in the determination of dg"^, by means of Eq.(70), can however 
allow only a qualitative estimate of triple parton correlations, while a better determination 
requires a dedicated study. Some preliminary results of these investigations were presented at 
the Conference "New trends in high-energy physics - Alushta(Crimea) - 2011" [63] 



6. CONCLUSIONS 

The aim of the present paper is to study the possibility of obtaining model independent 
information on multi-parton correlations, by measuring MPI in high energy collisions of hadrons 
with light nuclei. Two different kinds of correlations, in fractional momenta and in the transverse 
coordinates, are in fact unavoidably linked and cannot be disentangled, when studying MPI in pp 
collisions only. Already the simplest case, namely Double Parton Interactions in pD collisions, 
is characterized by novel and non trivial features, as compared to DPI in pp collisions. The 
component of the cross section, where both target nucleons contribute to the process, depends in 
fact only weakly on the correlations in transverse space. In addition, one has also a contribution 
from a interference term. All different contributions to the cross section can be disentangled 
and all the new unknown quantities appearing in the reaction and directly related to parton 
correlations, in fractional momenta and in the transverse coordinates, can be isolated in a way 
essentially model-independent. An interesting feature is that the interference term is expressed 
through the off diagonal parton distributions. By studying the interference term one may thus 
gain information also on the off diagonal parton distributions, in kinematical ranges not easily 
accessible through other processes. 

In order to disentangle triple parton correlations in fractional momenta and in the transverse 
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coordinates, one needs to measure Triple Parton Interactions in p or p ^He collisions. In 
fact, the knowledge of the cross section of Double and of Triple Parton Interactions in pp and in 
p D collisions is not sufficient to isolate all the unknown quantities, which appear in the reaction. 
On the other hand, keeping into account that the radii of D and of and ^He are not very 
different. Triple Parton Interactions in p D collisions can teach a lot on the contribution to Triple 
Parton Interactions in p or p ^He collisions, when only two target nucleons play an active 
role in the process. As discussed in Sub-Section 4.2, by taking advantage of the information 
on Triple Parton Interactions in pp and p D collisions, it's in fact rather simple to figure out 
how to obtain, from the cross section of Triple Parton Interactions in p or p ^He collisions, 
a model-independent, although only qualitative, indication on the different components of the 
cross section, with a direct link either with the triple parton correlations in fractional momenta 
or with the triple parton correlations in the transverse coordinates. 

Our conclusion is that MPI of hadrons with light nuclei have a great potential to provide model 
independent information on the multi-parton structure of the hadron. By measuring the cross 
sections with a given number of MPI on various nuclear targets, one may in fact identify different 
features of the incoming parton flux, which allows isolating diverse terms of the correlated multi- 
parton structure. To our knowledge, a result not accomplishable by other means. The option 
of studying MPI in collisions of protons with light nuclei at RHIC and to run, at some stage, 
light nuclear beams at the LHC could thus be highly rewarding, offering the possibility to exploit 
the remarkable potential of MPI in pA collisions to yield information on the many-body parton 
correlations and thus to provide unprecedented insight into the three dimensional structure of 
the hadron. 



Appendix A: The non-relativistic three-body wave function 

The nuclear systems we have considered {"^H, '^He) can be treated with a non relativistic 
dynamics in their center-of-momentum frame, but since they are involved in a highly relativistic 
process it is necessary to match this internal non relativistic dynamics with the overall relativistic 
treatment. To this end the original procedure used by Salpeter|64)|65j to reduce the Bethe- 
Salpeter equation to the Schrodinger equation will be followed as strictly as possible. We are 
able only to set a correspondence between nonrelativistic and relativistic wave functions, but 
not to build up a wholly deductive procedure as in the quoted refernces. For simplicity we treat 
both the constituent and the bound state as a spinless boson. The final aim is to use the non 
relativistic wave functions, as they are known from nuclear physics, in our relativistic calculation 
with the correct factors and the correct kinematical transformation. 

The starting point is given by an homogeneous equations in relativistic form as obtained by 
a Feynman graph representation, in term of the two-body scattering matrices as suggested by 
the Faddeev|66) |67| treatment of three-body t] scattering: 



U3{qi,q2,q3) = A(gi)A(g2) ^ / it^iqi + q2,k)Uj{qi - k,q2 + k,q3)dk, 
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+ G°3^3t3 Gl = A(gi)A(g2) 



(Al) 



together with the two analogous terms for Ui and U2- An iteration of the above equation shows 
that all the three terms Uj have in front the product of the free propagators of the constituent 
particles A(gi)A(g2)A(g3), as explicit in the graphical description in Fig. 10. 



u 



u. 



'2 U. 



Figure 10: Faddeev equation 



Defining T = qi + q2+ q-i the system is non-relativistic in the frame T = 0. 

When the two-body scattering matrices tj do not depend on the relative energies but only 
oh the three-momenta, as it happens in nonrelativistic dynamics, we can integrate Uj in ko- A 
generalization of this procedure suggests the Ansatz: 



Uj{qi,q2,qz) = «A(gi)A(^2)A(g3)$j(qi, q2, qs) 



(A2) 



This form has the meaning that, in the nonrelativistic limit, the bound particles are near the 
mass shell, so the singularities of the A— factors are the most important; morover, in this limit, 
antiparticles are not relevant, so the contribution of the antiparticle poles do not need to be 
taken into account in the integrations, explicitly: 



A(g) 



2oj Iq, 



1 



1 



■ UJ 



lo + 



1 



2uj qo 



■ UJ 



By inserting eq.(A2) into eq.(Al), dropping the common factors A((7i)A(q2) and integrating in 
ko we obtain: 



Afe) $3(qi,q2,q3) 



37 



= Afe) Yl / '^3(91 + Q2, k)A(gi - k)A{q2 + A:)$j(qi - k, qa + k, q3)dk (A3) 

-c-^ /■ 27r 1 o 

= Afe) > / t^^t;^^ ^ : t3(gi+g2,k)$j(qi-k,q2 + k,q3)d^/c 

^ J 4a;ia;2 (T - ^3)0 - ui - U2 + le 

The sub-energies are: = [(qi)^ + m^]^/^, lDi = [(qi — k)^ + m^]-^/^, 0)2 = [(q2 + k)^ + m^]^/^, 
We proceed by integrating over (53)0- Three sources of singularities need to be considered: 

- singularities of the propagators where ((73)0 appears directly; 

- singularities of the propagators containing {qi)o, (92)01 where (93)0 enters because we work at 
constant Tq. Actually it is convenient to make the position qi = (T — q3)/2 + l,q2 = {T — q3)/2 — l; 

- singularities of 1$, possibly originating from the two-body subsystem (1+2), which can be either 
poles like C/[(T — 53)0 — r/ + ie] (two-body bound states) or cuts like J p{ri)dr]/[{T — q3)o — ri + ie] 
(two-body scattering states). 

A direct inspection shows that the pole (^3)0 = 003 — ie has imaginary part with sign opposite to 
all other poles in (93)0, so its contribution gives the whole result of the integration. One obtains 



^3(qi,q2,q3) (A4) 

= 27rV / d^/c— Vt3(r-(?3,k)- 1 ^$^(q, _k,q2+k,q3) 

J 4wiW2 To - Ul - UJ2 - t03 + «e 

It is useful to make the positions: Tq = Mt = 3m + Eb, ujj = m + nj, where by construction 
it results Eb < 0, kj > . For the function defined as: 

¥'3(qi,q2,q3) = — ^3(qi,q2,q3 •) (A5) 

Eb - - K,2 - K3 

one obtains the equation: 

1 



2TTGoiEB) J ^3^jd^k with GoiEB) 



^3 

J^3 



Eb - Hi - K2 - K3 



which, together with the two analogous equations for ipi,ip2, represent the usual system of 
Faddeev equation for a non-relativistic three-body svstem|66]|67|. in fact Eb is the binding 
energy and the kj are the nonrelativistic kinetic energies. 

Defining as usual qi = T/3 -|- Ps/2 + 1, q2 = T/3 + ps/2 - 1, qg = T/3 - Ps, in the T = 
frame, the wave function depends on two three-vectors and setting: 



X]'^"'('ll'^2,q3) = ¥'(Ps,l) 



(f is the complete non-relativistic bound-state wave function, which can be made real and we 
must satisfy the normalization condition 



(f^d-'psd'H = 1 . (A6) 

The relations (A1,A4,A5) are linear and cannot give the normalization constant M. A nor- 
malization condition may be obtained by considering the total charge of the bound state. 
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Figure 11: Photon Tritium vertex 



The coupling of the Tritium with an external electromagnetic field is written as 

Jf'ik) = 2T^'f{k^) 

$(qi,q2,q3)A(gi)2(?fA(gi + A;) 

X A(g2)A(g3)$(qi + k, qa, q3)<5(T - ^ g) H 



(A7) 



In the limit A; — )• the zero component must give the total charge of the bound state j7o(0) = 
2To, or, in other words, the form factor must satisfy: /(O) = 1. Now we use the relations 
qi = T — q2 — qs and then, keeping the prescription of neglecting the antiparticle poles, we set: 



A{q)' 



1 



AqoUJ 



1 



qo - w 



The integration over the particle poles in j QSo^ performed in the frame T=0 with the result: 



joio) = -{2nfl n|7^(qi>q2,q3) 



1 



1 2 



To - UJl - UJ2 - LO-3, 



^(qi,q2,q3)5(J^q) =2Mt 

(AS) 

We compare (AS) with (A6), at first sight we might conclude that we need an M that depends 
on iOi, but we remember that in the present treatment we neglect k with respect to m, but not 
with respect to Eb, so we have finally 



TV = Imy^nMr/ 



IT 



(A9) 



The function ip has dimensions -3 in powers of q, see (A6), then the function $ has power zero in 
q as it must be, in fact $ plays finally the role of an effective coupling in a four-boson relativistic 
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vertex. Summing up the procedure: we start from ip as it is known from nuclear physics, we 
construct then eq.s (A5, A9) and from it we obtain, eq.(6), the expression of ^ that, after 
Fourier transformation, enters in the expression of the cross sections. 

The functions ip are given in terms of the three dimensional momenta in the cm. of the 
nucleus q^, we need them to be expressed in terms of the hght cone fractional momenta to 
evaluate the cross sections in the main text. In the case of interest two nucleons are on shell (see 
the figure). 




off shell 



on shell 



Figure 12: T/' 



^He vertex, one nucleoli is virtual and two nucleons are on shell 



Introducing the invariants 



Sl2 = {Pi +P2f 




(AlO) 



Sl3 = {Pl +Pzf 



{T-P2? 



one has 




2 / 2 




The light-cone 4-momenta components are given by 
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T = 



Pi = 



P2 = 



P3 = 



Ml 



, T_, 



m 



i± 

Zi(r_/3) 

,2 



, Zir_/3, pi_L 



m: 



2^ , Z2T_/3, P2±j 

, (3 - Zi - Z2)T_/3, -(pix + P2±) 



Z2(r_/3) 

+ (Pl-L + P2±)^ 

(3-Zi-Z2)(r_/3)^ 



(A12) 



where m?, = + p?, are the transverse masses. The 4-momentum conservation 



Mf = -^mi^^ + —mi^ + - — 



{P\ + (Pl± + P2±)]' 



(A13) 



implies 



p| - = (3 - Zi - Z2) 



Ml 
L~3~ 



i± 



m. 



2± 



3± 



Z2 3 — Zi — Z2 



(A14) 



In the T/^He centre of mass frame, the nucleon's energies E^, i = 1,2,3 are expressed in terms 
of the invariants U as follows 



El 
E2 
Ez 





+ - 


- S23 




2Mt 




M| 


+ - 


- S31 




2Mt 




Ml 


+ pI- 


S12 


2Mt 



(A15) 



The relations 



S23 

S31 



Mr. 



Mr. 



Z2 



Zi 
3 



3 '''^Z2 

allow to express Ei and E2 in terms of fractional momenta and transverse masses 



(A16) 



El 
E2 



1 



2-^1 , _2 



2Mt V ^ 3 ^-^^2/ 



(A17) 
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By using Ef = q| + m one obtains similar expressions for centre of mass 3-momenta qiz and 
q2z 



Keeping into account 



qiz 

q2z 



2Mt 
I 

2Mt 



3 

r2 ^2 



Mfi — ^ — m2_|_ 



^2' 



(A18) 



Mt = -El + £^2 + £^3, qi+q2 + q3 = 
one obtains the analogous relations for E^, and gs^: 



(A19) 
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2Mt L 



2 



m 



i± 



2± 
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Z^ 



Zx 



Zi 



m. 



2± 



(A20) 



Zi Z2 

which allow to express explicitly the non-relativistic nuclear wave functions as a function of the 
fractional momenta Zi and of the transverse momenta Pi±, i.e. in terms of variables invariant 
under longitudinal boost. 



Appendix B: Models worked out completely 

We present here two models where more explicite calculations are carried out after having 
introduced more or less strong simplifications of the real dynamics. 

1. A model with the Hulthen wave- function 



Since the Hulthen potential, with its relative wave function, is one of the simplest potential 
used in preliminary analyses of the Deuteron properties we present here a short derivation of 
some properties which are relevant for our investigations, in particular for an estimate of the 
relevance of the interference terms. 

The relative two body wave-function for the ground state (S-wave) , in r and in p representation 

is: 



T — K r '- 



Hp) 



7r(r — n) 



1 



1 



p2 + 



p2 + r2 



(Bl) 



where n = \fmE t = k, + the binding energy is £^ = 2m — Mo, If/i should represent the 
range of the potential, really it has been fitted phenomenologically to /x ~ 5k , m/2 is the reduced 
mass, assuming equal masses for the nucleons. The normalization is: 



4tt J \hir)\'^r^dr = 1 iir J \h{p)fp^dp = 
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Now we may calculate the mean value and the dispersion of the radial coordinate with the 
result: 

< ^ >= -7; — 7 ^ <r^ > - <r >^= p- — TTT^ (B2) 

2kt{k + t) [2tk{k + t)]^ ^ ' 

In the actual case the parameters satisfy the condition k, «t and there is a strong simplification: 

< r >Ri \/< > - < r >2 1/2k 
The longitudinal variable Z is studied in an analogous way and the results are: 



< Z^ >-- 



Mils 



<Z^ >-- 



16 



M4 L5 



16 



< p^ > +4m^ < p^ > 



(B3) 



with the mean values: 



< p^ >= KT < p^ >= KTiK^ + 3kT + r^) 

The qualitative behaviors of the parameters B, Z we learn are that the transverse extension is 
as large as one expected and the relative dispersion is quite large; the longitudinal variable Z is 
centered aroun 1, with a relatively small dispersion at least for k < r < m where we find: 

V< Z"^ > - < Z^ >7 < Z^ >= VSkt/Qtu 

The fact the the fractional momentum is slightly larger then 1 is due to the fsact the we are, 
really, considering an unsymmetrical situation where one of the bound nucleon is put on mass 
shell, the configuration is symmetrical in the space variables p, but it is not symmetrical in the 
relative energies. It seems that the more interesting result is the dispersion in Z, in fact this is 
the parameter which says how much the contribution of the interference term, where Z 7^ Z', 
differs from the diagonal term. 

We are also interested in a mixed representation where the transverse degrees of freedom 
are expressed in space variables B while the longitudinal degree is given in light-cone variables 
p+,P-- We recall that we are interested in a particular kinematical situation where one of the 
bound nucleon is treated as real on mass shell but we still are on the center-of-momentum frame 
so that the two three momenta are opposite, then we shall consider a longitudinal boost. In this 
situation it results pz = m]_/2p- — p_/2 = + p^ With these definitions we obtain: 



where 



1 



P- 



p^ + 



1 



1 



p'j_ + w"^ 



p\ + v"^ _ 



p 



L_ + 2p_ \/ — + = — 2p_ ^/ — + rr? 



and OL is either k or r. The definition Z = 2p-/D- gives in the center-of-momentum frame 
P- = M£)Z/2 and so we get: 



h{B,Z) = —^kt{k + t) 



-\Ko{w^B) - Ko(v^B) 
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Figure 13: Graph of h{B = 0) as function of p^/m. Arbitrary vertical normalization; the upper curve 
corresponds to fi = 3k, the lower curve corresponds to /i = 5k. 



Z 



= [Ko{WrB) - Ko{VrB)] 



(B4) 




In order to apply these expressions to the Deuteron's case we take into account that n < fi < 
m; we have already noted that we are interested in small values of B compared with the nuclear 
scale 1/2k, so we look for the limit -B — )• which gives Ko{wkB) — Ko{vk.B) — >• \q.{v^/wh). 

We can study numerically the above limiting form of /i, the result for E/m = 0.0021 are 
presented in the graph as function of u = p-/m., which is slightly different from Z = 2mu/Mo- 
Since ^ = 5k is a phenomenological fit without direct interpretation the numerical study has been 
performed also for fj, = 3k which would better describe a binding potential generated by pion 
exchange; the qualitative conclusions are however the same with both choices of the parameters. 
In Fig. 13 we plot h{0, Z) as a function of Z in the cases fi = 3k and fi = 5k. A numerical study 
shows that the shape of h{B, Z) has little variation with increasing B. In accordance with the 
previous results the spread in Z is sizable, in order to go from the maximum of h to half of this 
maximum Z must vary at least by 0.1. So the product of two function with Z different of some 
per cent, and then its integral, is not very depressed with respect to a diagonal term. This in 
practice means that the interference terms are smaller than the diagonal ones, but not by orders 
of magnitude. 



Here the simplification is heavier, in the aim of dealing explicitly with the transverse degrees 
of freedom for all the graphs we considered. No correlations among transverse and longitudinal 
degrees of freedom are taken into account, moreover the transverse distributions are Gaussian. 
One of the results, which seems us of more general validity is that the transverse dynamics is. 



2. 



A Gaussian model for the transverse dynamics 



a. General features of the model 
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in comparison with the longitudinal, less sensible to the difference between diagonal terms and 
interference terms. 

The one-parton inclusive distribution is written as: 

Ti = G{x)h{b) h{b) = ^ expi-b^/R"] . (B5) 
The two-parton distribution is: 



T2 = KG{xi)Gix2)f2{bi,b2) h{biM) = ^^^2)2(i _ ;^2) exp[-(6^ + bl + 2\bib2)/R\l - X") 



I 



h{b^M)db2 = h{hi) (B6) 



K controls the parton multiplicity, A their spatial correlation; both can still depend on the 
fractional momenta Xj. 

The three-parton distribution, with the minimal number of new parameters, is written as: 

Ts = G{xi)G{x2)G{x^)KMhiMM) 

1 r (i-A)E.^^ + 2AE.<,(6,-6,) . 

^'^^'^^'^ = {.R-f{l-2X){l + Xf [ i?2(l + A)(l-2A) J ^^^^ 

where we have implemented the requirement that /a be symmetric in the parton variables 
and that 



/ 



h{bi,b2M)db2, = f2{bi,b2) ■ 



The information embodied in the parameters and A can be seen in this way. The distribu- 
tions /at are inclusive, if the corresponding exclusive distributions, integrated on the longitudinal 
variables in a given interval, were Poissonian, then Kj^ would be 1. Of course even a Poissonian 
distribution of the partons' multiplicities can imply spatial corralations if A 7^ 0. 

The nuclear distribution, i.e. the square of the wave function are expressed for the deuteron 

as: 

g{z- b) = W{Z)F{B) F{B) = exp[-SV^^] 0<Z<2 B = Bi - B2 (B8) 
Since the distribution will be at the end integrated, for the Tritium we write: 



g{Zi, Bi) = W{Zi, Z2, Zs)S{Zi + Z2 + Zs)F{Bi) 0<Z<3 

^(^i) = [ - - B2f + {B2 - Bsf + {B2 - BsfysS^] (B9) 

The normalization is 

J W{Z)dZ =1, j W{Zi, Z2, Z^)5{Zi + Z2 + Z3)dZidZ2dZ3 = 1 
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In nonrelativistic case Z The transverse three-body distribution satifies 

^ F{Bi)dBs = expl-B^S^] B = Bi - B2 , 

which defines the normaUzation, it has been aheady noted that the sizes S in the deuteron and 
in the Tritium, although similar, are in fact different. 



b. Free nucleons 



As reference quantities we consider, within the model, the simple double and triple hard 
scattering among free nucleons, at fixed fractional momenta. The simple hard scattering is 
described by 

= a^^.G{x)G{x') J h{b)h{b - (3)dbdp = a^^,Gix)G{x') 
due to the normalization of the transverse distributions. The double hard scattering is described 

by 



a2{xi,X2;x[,X2) = Klax^^>_G{xi)G{x'-i)a^^^>^G{x2)G{x'2)Io 

lo = j f2ibi, 62)/2(6i -P,b2- P)dbidb2dp (BIO) 

The triple hard scattering is described by 

(73 {xi,X2,X3;x[,X2,x's) = Kla^^^>^G{xi)G{x\)a^^^>^G{x2)G{x'2)cF^^^>^^^ 
Jo = j Mh, b2, 63)/3(6i -P,b2- 63 - mbidb2db3dp (Bll) 

All the transverse variables "5" are two-dimensional and the cylindrical symmetry is always 
preserved, so all the calculations regarding the transverse variables have this standard form 



G 



d y exp[-y ■ M ■ y] 



= C 



detM 



(B12) 



Here M is a iV x iV matrix, C embodies the normalizing factors as given in Eq.s (B5, B6, B7) 

In the cases considered above, after a rescaling of the variables, b bR-\/l — in the first 
case and b bR\/l + A\/l — 2A in the scecond case the matrices take the form 



Ml 



2 

2A 
-1- A 



2A -1-A 
2 -1-A 
-1-A 2 + 2A 



(B13) 



Ml 



2(1 -A) 2A 2A -1 

2A 2(1 - A) 2A -1 
2A 2A 2(1 -A) -1 



(B14) 



-1 



-1 3(1 + A) 
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So that we have, at the end: 



1 



1 



Jo 



47ri?2 1 + A 
1 



12(7ri?2)2 1 + 4A + 2A2 



(B15) 



c. Bound nucleons 



We start considering the scattering process when one of the nucleon is bound in a deuteron 

aD,i{x,x') = a^^:G{x)G{x' /Z)W{Z)dZdxdx' j h{b)h{b - Bi)F{Bi - B2)dbdBi 

By integrating in B2 the factor F one gets simply 1. According with the discussion in Section 
3, in the distribution W the variable Z is shrunk around the value Z = \ so that we take 
J" G{x' lZ)dZW{Z) Ri G{x') J dZW{Z) = G{x') and we obtain the same expression as for the 
free case. Prom the simple procedure described is is seen that the same result holds for the case 
of the Tritium, with a suitable redefinition of the size S also when we consider double or triple 
hard scatterings. 

When we look at hard scatterings where more than one bound nucleon participates new 
features appear. In the simplified model here described the longitudinal degees of freedom are 
factorized we have anyhow seen that Z ^ 1 approximation is inconsistent with the conservation 
of longitudinal momentum in the nondiagonal cases, but we do not have anything new to add to 
this point: we investigate the transverse degrees of freedom of the partons (inside the nucleon) 
and of the nucleon (inside the nucleus) which arc dynamically connected. 

Now we study the double scattering when both nucleons of the deuteron are involved: (here 
and below the x, x' arguments of a will be usually omitted) 

It has been shown that in this case there are two possibilities: the direct term and an inter- 
ference term. In this last case we cannot work simply with the density distribution of partons 
but we need the "wave function", whose absolute square gives the distribution of the partons 
inside the hadron. In our case, where the distributions are Gaussian we take as "wave function" 
simply the square root of the distribution, i.e. we ignore the possibile phases. Then the relevant 
quantities are: 



ll = J f2{biM)h{bi - Bi)h{b2 - B2)F{Bi - B2)dbdB 

li = j f2{biMW h{hi - Bi)^h{bi - B2Wh{b2 - B2Wh{b2 - Bi)F{Bi - B2)dbdBBlQ) 

The subsequent calculations are very similar to the previous one. There is only a new feature, 
the new dimensional parameter S. We have an integration over the transverse variables b and 
we define adimensional variables through b — >■ bR and, for convenience we introduce v = B? / S'^. 
The Gaussian integration involves the calculation of 4 x 4 determinants and finally, taking into 
account the normalizations we get, for the diagonal and for the interference terms, respectively: 

" 7r[52 + 2(2 + X)B?] ^2 = ;r(2 + A)[52 + 2i?2] ^^^^^ 
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In the double collision involving a Tritium there is necessarily a spectator. Since by integrating 
the three-body distribution (B5) over the spectator's coordinates we obtain the two-body distri- 
bution (B4) the final expression is the same, provided we rescale, in (Bll) the size according to 
the experimental values. 

Now we consider that triple hard interaction where one of the bound nucleons interacts twice, 
another only once and there is, in the Tritium case, a spectator. It has already been shown that 
there are two kind of processes in this case, a diagonal term and some interference terms which 
must be treated separately. 

We do not repeat the consideration about the longitudinal variables, for what concerns the 
transverse variables the integration implies the calculation of 5 x 5 determinants and the final 
result is, for the diagonal term: 

'^^ ^ T+A 47r2i?2[S2 + (3 + A)i?2] ^^^^^ 



For the interference terms we get 



T»,2 ^ 

2 (1 + A)(6 - A) 7r2i?2[(4 - A)S2 + 2(4 - 2A - \^)B?\ 



A 



■^2 



1 - A 1 



(1 + A)(3 + A) 7r2E2[(3 _ A)52 + 4(1 - X)E?] 

jiA ^ 1 ~ ^ ^ Co-,Q^ 

2 (l + A)(4-2A-A2)7r2i?2[752 + 2(6-A)i?2] ^'^^''^ 

Finally we consider the process where we have three bound nucleons, all interacting once, evi- 
dently now we must consider a Tritium. We find a diagonal term and two different interference 
terms, the integration implies the calculation of 6 x 6 detrminants and the final result is, for the 
diagonal term: 

" 3 7r2[52 + 2(2 + A)E2]2 (^^^^ 
For the two interference terms we get: 



ji,2 



1 



3 7r2(2 + A) [(52 + 2i?2)(52 + 2(2 + X)R'^] 

Ti,3 ^ 4 16 

^ 3 7r2[(7 + 3A)S2 + 8(2 + A)i?2]2 ^""^^^ 



A feature that appears very clearly from the model but that reflects a more general property 
is the dependence on the geometrical parameters: since the cross section has dimension i"^ and 
it embodies three factors a, with dimension £2 each one there must be a factor When only 
one nucleon interacts this factor is necessarily 1/i?^, when two nucleon interacts the factor is 
l/i?^(/xS'^ -I- vR'^), when three nucleons interact the factor has the form l/(/i-S'^ + vR'^)'^. The 
meaning appear clear considering the hypothetical situation S » R. In this case the first cross 
section remains unaltered, the second vanishes as l/S'^, the third one vanishes as 



Appendix C: Infrared behavior 



We give a short look at the infrared properties of the amplitudes and densities we have used. 
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Using eq.s (5,9) we would conclude that the one particle density has the behavior T{z; b) oc 
1/z, but we know that this is not supported by experimental data, so we must conclude that 
the vertex (p must show also an infrared singularity in such a way that r{z;b) oc l/z^+''. Since 
the integration (with an infrared cutoff) of the two-body distribution must give the one-body 
distribution the same behavior must be found in the two-body vertex cj). 

Thus we find a relevant simplification of the two-parton amplitude when one of the two 
partons has its four momentum very soft, it is the particular case of the general features found 
in the emission of soft particles. For dcfinitcncss we consider the free proton case: when the 
parton four-momentum goes to zero linearly in all its components, we have li± oc xi. Then in 
the expression for ip2 we neglect the term in and also xi with respect to X2 with the result. 

, i 

V2L+XI X2[m? - Mil {I - X2)] - hi 

but we have seen that ^ must have a singularity, we extract it by writing (f) oc In the 

same kinematical configurations also the other factor used in defining the densities V can be 
decomposed as 

Xi,X2 X2 

~Xl- Xl«X2. 

I — Xl - X2 I — X2 

So the original density is decomposed into two factors; the second factor is precisely the proba- 
biliy of finding one parton with finite fractional momentum, as seen in eq.(5),the first one can be 
thought as the usual infrared term of QED correced phenomenologically enhancing the singular- 
ity; the well known term of QED is P ■ e/P ■ q, it could represent also a gluon emission, but since 
we have not taken into account the spin we do not have the numerator P ■ e; the denominator is 
dominated by the "large" component xiL-^. 

We use the mixed representation of tp with longitudinal momenta and transverse coordinate, 
so we should perform the Fourier transformation in l±. The operation on l2± yields precisely 
ip{x2; 62)) on Zij^ the integration must run only on the infrared domain of Zix, limited by a cut off 
|^i±| < ii.r. for dimensional reasons the result is proportional to a function G(6i4.r.) and we 
get for the density r(xi, X2; 62) — •^i,,..G'(&i£i.r.)r(a;2; 62) /2.ti. Integrating over the longitudinal 
infrared momentum xi we obtain the usual phenomenological divergence dx/x^'^'^. 

It is in fact a very general property that the soft emission is relatively independent by the rest 
of the dynamics, so it holds also for the more complicated expressions, like the terms Wj and also 
for three-body amplitudes, which become completely factorized in the limit xi « X3; X2 « 
X3. In detail the treatment applies directly to a soft emission originating directly from a free 
nucleon. When the parton emission comes from a nucleus, one finds the usual complication due 
to the binding, but when the treatment is applied to a configuration with z « 1, the previous 
discussion still holds, see eq.21. The factorization which is produced in this particular case can 
be useful in order to make a simpler relation, at least in this kinematical configuration, between 
the diagonal and the interference term. 
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